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Abstract 

A  commonly  accepted  airborne  phased  array  radar  model  simplifies  the  ana¬ 
lytical  derivation  by  assuming  a  waveform  is  perfectly  matched  in  range  and  Doppler 
shift.  This  assumption  means  the  matched  filter  output  is  effectively  constant  for 
all  possible  received  scatterer  Doppler  and  range  mismatches,  greatly  simplifying 
the  analytical  development  from  that  point  forward.  This  research  removes  the 
matched  Doppler  and  range  assumption  and  examines  the  effects  of  several  common 
waveforms  on  the  model’s  fidelity  along  with  the  associated  impact  on  radar  perfor¬ 
mance,  both  non-adaptive  and  adaptive.  Analysis  is  completed  using  power  spectral 
density  comparisons  and  the  fully  adaptive  output  signal  to  interference  plus  noise 
ratio  comparison.  Results  indicate  that  the  model’s  fidelity  is  impacted  little  by  the 
Time  Frequency  Auto  Correlation  Function,  ffowever,  change  in  bandwidth  from 
the  compressed  waveforms  does  impact  the  model.  Increased  bandwidth  introduces 
more  thermal  noise  which  dominates  clutter  returns.  Therefore,  the  clutter  problem 
becomes  less  difficult.  The  trade-off  is  a  reduction  in  the  resolution  capability  of  the 
clutter  spectrum. 
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An  Examination  of 


Range  and  Doppler  Mismatch 
AND  THEIR  EFFECTS  ON  RADAR  MODELING 


I.  Introduction 

Radar  models  and  simulations  provide  insight  into  the  operation  of  radar  in  real 
world  conditions.  A  robust  simulation  is  created  by  generating  a  complete 
radar  model.  The  complete  radar  model  includes  accurately  representing  the  physical 
radar  with  realistic  operating  parameters  and  environment  characteristics. 

Assumptions  are  a  tool  used  in  many  models  and  simulations  to  simplify  diffi¬ 
cult  problems.  When  creating  the  basic  radar  model,  outlined  later,  assumptions  are 
made  to  simplify  the  mathematics  and  simulations.  The  majority  of  these  assump¬ 
tions  are  valid  under  normal  operating  circumstances  and  have  no  adverse  effects  on 
simulation  results  when  compared  to  real  world  operations. 

There  are  two  assumptions  in  focus  throughout  this  research.  The  assumptions 
are  that  there  is  no  mismatch  in  range  and  that  a  Doppler  tolerant  waveform  is 
transmitted,  even  if  the  actual  waveform  is  not  Doppler  tolerant.  This  research 
explores  the  removal  of  these  assumptions  and  analyzes  the  effects  of  Doppler  and 
range  mismatch  on  the  radar  model. 

1.1  Background 

The  first  step  to  examining  this  issue  is  creating  a  radar  model,  commonly  ac¬ 
cepted  by  many  in  the  Space  Time  Adaptive  Processing  community,  for  comparison. 
This  model  incorporates  a  Time  Frequency  Auto  Correlation  Function  (TFACF)  of 
the  waveform,  a  formed  antenna  beam,  beam  steering  abilities,  and  any  effects  from 
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thermal  noise  as  well  as  clutter  returns.  The  TFACF  is  where  the  assumptions  are 
made. 

The  next  step  is  to  mathematically  derive  the  TFACF  expression  for  the  non- 
Doppler  tolerant  waveform.  The  TFACF  takes  into  account  possible  shifts  within 
the  matched  filter  due  to  Doppler  and  range  cell  mismatch.  A  limitation  of  no  more 
than  plus  or  minus  a  half  chip  width  shift  is  placed  on  the  derivation.  This  limitation 
allows  the  target  Signal-to-Noise  Ratio  (SNR)  to  remain  constant.  A  new  weighted 
TFACF  expression  fits  into  the  regular  radar  model  at  the  same  location  as  the 
Doppler  tolerant  expression  of  the  TFACF. 

Using  the  TFACF  expression,  simulations  are  run  for  several  common  wave¬ 
forms.  These  waveforms  are  simulated  using  both  radar  models;  the  commonly 
accepted  radar  model  (Space  Time  Adaptive  Processing  community)  which  assumes 
Doppler  tolerance  and  perfect  range  matching,  and  the  modified  radar  model  which 
implements  the  TFACF  allowing  for  Doppler  and  range  mismatching.  Results  from 
each  model  are  compared  against  each  other  using  Power  Spectral  Density  (PSD) 
methods  as  well  as  the  fully  adaptive  radar  processor  Signal  to  Interference  plus 
Noise  Ratio  (SINR)  method.  Analyzing  the  results,  conclusions  are  drawn  as  to 
whether  or  not  these  effects  are  necessary  for  future  modeling. 

1.2  Organization 

In  chapter  two,  the  commonly  accepted  radar  model  is  outlined  showing  the 
development  of  the  model  including  aperture  geometries,  transmitting/receiving  sig¬ 
nals,  steering  vector  roles,  and  the  environment  characteristics.  Within  this  devel¬ 
opment,  the  assumptions  under  question  are  highlighted.  Then,  in  chapter  three, 
the  TFACF  expression  that  accounts  for  Doppler  and  range  shifting  is  found  as  well 
as  how  they  relate  with  in  the  radar  model.  Chapter  four  shows  the  results  and 
provides  insight  and  analysis  on  how  the  two  radar  models  differ.  A  summary  of  the 
entire  research  is  covered  in  chapter  five. 
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II.  Literature  Review 


A  complete  understanding  of  the  radar  model  is  needed  before  addressing  any 
effects  that  Doppler  and  range  mismatching  creates.  This  chapter  outlines  each 
separate  segment  that  encompasses  the  radar  model,  along  with  signal  processing 
techniques  and  overall  system  assumptions. 

The  radar  model  has  been  developed  and  modified  throughout  the  years.  This 
development  of  the  radar  model  presented  in  this  chapter  is  based  on  previous 
work  [1,  5,  7]  and  written  in  this  authors  own  words.  The  radar  being  simulated 
is  based  on  an  airborne  phased  array,  pulsed  Doppler  system  illuminating  a  tar¬ 
get  point  scatterer.  Therefore,  all  characteristics  and  properties  of  a  pulse  Doppler 
system  are  accounted  for  in  the  model. 

Combining  many  different  aspects  together  creates  the  radar  model.  The  first 
aspect  is  to  understand  the  geometry  of  the  aperture.  The  geometry  effects  the 
antenna  transmit  pattern  which  is  necessary  for  the  radar  range  equation.  The  next 
step  is  to  examine  the  transmit  signal.  After  the  transmit  signal  is  reflected  from  a 
point  scatter,  the  received  signal  is  examined  and  processing  occurs.  A  by-product 
of  this  processing  is  the  steering  vector.  This  steering  vector  allows  for  weighting 
on  the  antenna  beam  to  allow  for  steering  the  beam  in  space  and  time.  The  next 
section  of  the  model  is  builds  the  clutter  and  noise  environment.  However,  before 
building  the  noise  and  clutter  section  of  the  model,  the  element  and  array  pattern 
is  needed.  These  patterns  are  needed  because  the  clutter  model  is  dependant  on 
the  radar  equation  which  includes  the  antenna  array  and  element  patterns.  At  this 
point,  the  model  is  almost  complete.  The  last  two  sections  of  the  radar  model  include 
interferers  with  the  received  signal.  Thermal  noise  occurs  within  each  element  in  the 
array  and  needs  to  be  accounted  for  and  modeled.  The  last  interferer  is  ground 
clutter. 
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A  Azimuth  Elements 

Figure  2.1:  This  shows  the  array  configuration.  Antenna  ele¬ 
ments  are  placed  in  the  x  and  z  direction  with  element  spacing 
of  dx  and  dz.  Element  A"*  is  located  at  (x,z)  =  (0,0). 


2.1  Aperture  Geometries 

The  aperture  structure  is  the  first  block  to  be  examined  in  the  model.  Under¬ 
standing  the  geometry  is  necessary  for  development  later  in  the  model.  Element  and 
array  patterns  along  with  the  steering  vector  are  directly  dependant  on  the  set-up 
of  the  aperture  geometry. 

The  antenna  array  structure  consists  of  elements  with  the  capacity  of  trans¬ 
mitting  and  receiving.  Elements  placed  along  the  x-axis  as  well  as  the  z-axis  create  a 
rectangular  array.  Each  element  is  uniformly  spaced  at  intervals  of  dx  and  dz.  How¬ 
ever,  dx  and  dz  are  not  necessarily  equal.  Notational  wise,  there  are  N  dx  elements 
in  the  ^-direction  (azimuth  elements)  and  P  represents  elements  in  the  ^-direction 
(elevation  elements).  Figure  2.1  shows  the  geometry  of  the  array. 

A  generic  expression  describing  the  location  of  the  npth  element  is  given  by 


dnp  =  —ndxx  -  pdzz , 


(2.1) 
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where  n  ranges  from  0  to  N  —  1,  p  ranges  from  0  to  P  —  1,  and  x  and  5  are  unit 
vectors  along  the  x-and  z-axis,  respectively. 

The  radar  coordinate  system  uses  9  to  describe  elevation  and  (p  as  azimuth. 
The  reference  element  is  the  upper  left  of  the  array  (annotated  in  Fig  2.1  by  A"*). 
Thus,  9  is  negative  when  measured  from  bore  sight  toward  the  ground.  Also,  by 
placing  the  reference  element  at  this  position,  dnp  becomes  negative  in  both  the 
x  and  ^-direction.  Conversion  from  the  radar  coordinate  system  to  the  Cartesian 
coordinate  system  (needed  later  for  the  antenna  array)  is  given  by 

k(9,  (p)  =  cos  6  sin  <px  +  cos  9  cos  (py  +  sin  9z.  (2.2) 

The  unit  vector  k  is  a  reference  vector  in  the  direction  of  the  target  described  by  9 
and  (p ,  yet  expressed  in  Cartesian  coordinates. 

Establishing  the  aperture  geometry  allows  for  continuing  development  with 
the  rest  of  the  model.  In  order  for  a  radar  to  be  able  to  detect  targets  in  space,  the 
radar  must  first  transmit  a  signal.  Therefore,  an  examination  of  the  transmit  signal 
is  needed. 

2.2  Transmit  Signal 

The  radar  transmit  signal,  labeled  ‘A’  in  Fig.  2.2,  is  a  sinusoidal  waveform, 
mathematically  expressed  as 


s(t)  =  atu{t)e^ot+^ ,  (2.3) 

where  at  is  the  transmitted  pulse  amplitude  and  u(t)  is  the  envelope.  The  com¬ 
plex  exponential  describes  the  carrier  frequency  u)0  and  the  random  phase  shift  ip. 
The  random  phase  shift  is  necessary  since  the  phase  state  of  the  transmitted  cosine 
waveform  is  not  known  at  the  time  of  transmission. 
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Figure  2.2:  This  figure  represents  the  overall  setup  of  the 

radar  model. 


Decomposing  u(t)  provides  more  information  about  the  pulse  train.  This  de¬ 
composition  is  given  by 

M—l 

”  (2.4) 


M—l 

Mm (£)  =  ^  u{t  —  mTr ), 


m= 0 


where  the  summation  of  individual  pulses  u{t)  comprise  the  complete  envelope  Um{E) 
and  Tr  is  the  Pulse  Repetition  Interval. 


An  examination  of  the  energy  in  the  transmit  signal  provides  an  important 
check  later  in  the  development.  Knowing  the  amount  of  energy  in  the  transmit 
signal  creates  a  ceiling  for  the  received  signal.  If  the  energy  in  the  received  signal  is 
greater  than  the  transmitted  signal,  then  there  is  an  error  within  the  mathematical 
calculations  or  there  is  an  outside  source  providing  energy  to  the  system  to  spoof 
the  radar.  However,  for  this  model,  no  spoofing  is  added  to  the  system.  Therefore, 
checking  the  energy  in  the  returned  signal  is  a  good  error  check  for  the  mathematics. 


The  energy  in  the  transmit  signal,  normalized  over  one  Pulse  Repetition  Inter¬ 
val  (PRI),  is  defined  as 


fmTr 


Em  = 


\s(t)\2dt, 


(2.5) 


6 


where  the  subscript  M  indicates  the  number  of  pulses  comprising  the  total  energy 
and  s(t)  is  the  transmitted  signal.  Allowing  Em  to  be  defined  only  over  a  single 
pulse  (M  =  1),  the  fundamental  energy  becomes 


T 


(2.6) 


where  T  is  the  pulse  width.  Equation  (2.6)  is  then  normalized  to  achieve  an  energy 
value  of  one  shown  by 

u(t)\2dt  =  l.  (2.7) 

Now  adding  the  amplitude  at,  the  energy  in  a  single  pulse  becomes 

E  =  ^  [  \atu(t)\2dt.  (2.8) 

1  Jo 

Since  at  is  not  a  function  of  t,  a 2  is  pulled  out  of  the  integral.  This  action  simplifies 
the  energy  equation  for  a  single  pulse  to 

E  = [r  \u(t)\2dt  =  al  (2.9) 


Equation  (2.7)  states  that  this  integral  equals  one.  Therefore,  the  energy  in  a  single 
pulse  with  amplitude  at  must  equal  a2.  Then,  due  to  integration,  the  total  energy 
in  the  transmitted  signal  is  scaled  by  a  factor  of  M,  the  number  of  pulses.  Modify¬ 
ing  (2.5),  the  total  energy  in  M  pulses  is 

/»MTr 

Em  =  /  \s(t)\2dt  =  Ma2.  (2-10) 

Jo 


The  transmit  signal  is  now  defined  and  some  energy  constraints  are  placed 
on  the  signal  for  error  detection  later  in  the  development.  A  point  scatter  reflects 
the  transmit  signal  back  in  the  direction  of  the  radar.  This  returned  signal  carries 
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information  about  the  point  scatterer.  To  gather  this  information,  the  received  signal 
needs  examination. 

2.3  Received  Signal 

The  received  signal,  labeled  ‘B’  in  Fig.  2.2,  is  the  same  as  the  transmitted 
signal  with  two  differences.  The  first  difference  is  an  added  Doppler  shift  /d  due  to 
the  relative  velocity  between  the  aircraft  and  the  target.  The  second  difference  is 
a  time  delay  to  the  npth  element,  rnp.  Incorporating  these  differences,  the  received 
signal  r(t)  becomes 

rnp(t)  =  s{t  -  Tnp)ej2^-T^.  (2.11) 

An  assumption  is  made  here  that  all  of  the  antenna  elements  are  in  the  far  field 
and  therefore  the  Doppler  shift  variation  across  the  elements  is  negligible.  Assuming 
negligible  Doppler  shift  simplifies  the  math  later  in  the  development.  The  time  delay 
is  needed  due  to  the  returned  signal  striking  the  array  in  a  non-parallel  fashion. 
Tracking  time  delays  is  necessary  to  be  able  to  coherently  add  the  phases  and  steer 
the  array. 

Before  proceeding  to  the  development  of  the  received  signal,  there  is  a  condition 
that  must  be  considered.  This  condition  is  Doppler  foldover  or  Doppler  ambiguity. 
Doppler  foldover  occurs  when  the  Doppler  shift  is  greater  than  the  Pulse  Repetition 
Frequency  (PRF)  fr.  Therefore,  as  a  check  for  ambiguity,  the  Doppler  frequency  is 
normalized  by  the  PRF, 


where  u  represents  the  normalized  Doppler  frequency.  If  u;  >  || | ,  then  there  is 
Doppler  ambiguity.  Otherwise,  the  Doppler  frequency  is  unambiguous. 

Looking  back  at  the  received  signal  in  (2.11),  there  is  a  time  delay  compo¬ 
nent  that  requires  understanding.  To  find  the  total  signal  delay  to  each  element,  a 
summation  of  the  total  round  trip  time  of  the  signal  (Tr)  and  the  delay  from  the 


reference  element  to  the  npth  element  (r  )  is  performed, 


Tnp  /  / 1' 


T, 


np ' 


(2.13) 


The  reference  time  delay  requires  further  attention.  t  describes  the  relationship 
between  the  npth  element  to  that  of  the  reference  element.  This  distance  is  otherwise 
known  as  the  differential  delay  [1],  This  delay  is  a  projection  of  the  unit  vector  in 
the  direction  of  the  target  (described  in  terms  of  6  and  (ft)  onto  the  vector  describing 
the  location  of  the  element  from  the  reference  element  (dnp)  divided  by  the  speed  of 
light  c.  This  projection  results  in 


T  = 
'  np 


k(9 ,  (ft)  ■  d. 


np 


(2.14) 


Expanding  k  and  dnp:  rnp  becomes 

/  [cos  6  sin  <ftx  +  cos  9  cos  (fty  +  sin  9z\  ■  [— ndxx  —  pdzz] 

T  —  - 

np  „ 


(2.15) 


The  array  exists  only  in  the  x  and  z-axis,  therefore  the  y  term  goes  to  zero  resulting 
in 

/  —  ndx  cos  6  sin  (ft  —  pd-,  sin  (ft 

TnP  = - - - - - •  l2-^) 

Now  that  the  time  delay  is  known,  a  conversion  of  the  time  delay  into  the  spatial 
domain  is  necessary  due  to  a  phase  delay  acting  on  the  original  and  Doppler  shift 
frequencies.  To  convert  (2.16)  into  the  spatial  domain,  both  sides  of  the  equation 
are  multiplied  by  u0 

,  —ndx  cos  6  sin  (ft  —  pdz  sin  (ft 

uj0rnp  =  uj0 - .  (2.17) 


9 


To  help  simplify  (2.17),  the  variables  dx  and  dz  and  are  assigned 

dx  cos  9  sin  <p 


and 

dz  sin  6 

By  assigning  the  variables,  (2.17)  becomes 

UoT^p  =  -2tt  (mfy  +  pdz).  (2.20) 

With  a  better  understanding  of  t  and  issues  regarding  Doppler  foldover, 
development  of  the  received  signal  in  (2.11)  continues.  The  transmitted  signal  in  (2.3) 
is  inserted  into  (2.11)  producing 

rnp(t)  =  aruM(t  -  Tnp) e-7271^0 9~Tnp) e^fdit-Tnp) ejip ,  (2.21) 

where  ip  is  the  random  starting  phase,  ar  is  the  received  signal  amplitude,  and  fQ 
is  the  transmit  frequency.  The  parameter  ar  will  be  much  smaller  than  at.  This 
reduction  in  signal  strength  is  due  to  radar  range  equation  effects.  Substituting  the 
total  time  delay  in  gives 


(2.18) 

(2.19) 


rn,{t)  =  a,.uM(t  -Tr-  T^)e^<'"+«<‘-T*-Oe*‘.  (2.22) 

At  this  point  another  assumption  is  made.  This  radar  model  is  only  for  narrowband 
operations.  This  assumption  means  the  time  delay  Tnp  to  each  element  is  much 
smaller  than  the  pulse  width.  Therefore,  the  time  delay  within  the  unit  pulse  of  the 
signal  can  be  ignored.  If  the  transmitted  pulse  was  a  wideband  pulse,  there  would 
be  decorrelation  across  the  face  of  the  array  which  would  need  to  be  accounted  for. 
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This  assumption  reduces  (2.22)  to, 


rnp(t)  =  aruM(t  -  TR)ej2^fo+fd)it-TR-T'^e^ 


(2.23) 


Expanding  the  exponential  terms  produces 


f  np 


(t)  =  aruM(t  -  rfi)ej27r(/°+/d)te~i2^(/°+/d)Tflp~J^(/°+/d)T^ 


(2.24) 


Examining  the  bracketed  term,  there  is  a  <l'Kf0Tnv  term  which  is  also  written 
as  uJoTnp .  Substituting  (2.20)  in,  the  bracketed  exponential  term  is  written  as 


ej2-KfdTnp  ej  2rr(m}x  +ptfz ) 


(2.25) 


Rewriting  (2.24)  with  this  expanded  exponential,  the  received  signal  becomes 


At  this  point,  all  the  bracketed  exponential  terms  in  (2.26)  are  not  a  function  of 
time.  Therefore,  these  terms  maintain  a  constant  value  over  the  duration  of  the 
pulse  and  are  rolled  into  the  ijj  term.  No  information  is  lost  in  doing  combing  these 


terms.  Later  in  the  development,  the  e ^  term  and  all  the  other  exponentials  that  it 
encompasses,  is  taken  into  account.  The  received  signal  now  equals 


(t)  =  aruM{t  -  TR)e^{fo+Sd)tej2^n^+pi)z)e^. 


(2.27) 


At  this  point  in  the  model,  the  signal  in  (2.27)  has  been  received  at  the  element 
and  is  ready  for  signal  processing.  There  are  steps  that  happen  at  this  point  in  real 
world  operations  that  are  not  simulated.  First,  the  received  signal  is  passed  through 
a  low  noise  power  amplifier.  Then  the  signal  is  sent  through  a  coherent  oscillator, 
labeled  ‘C’  in  Fig.  2.2.  The  oscillator  strips  the  carrier  frequency  from  the  received 
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signal  leaving  the  signal  at  baseband.  After  these  two  steps,  the  received  signal  now 
looks  like, 

rnp{t)  =  aruM(t  -  TR)ej2nfdtej2n^+^^e^.  (2.28) 

The  next  step  in  the  signal  processing  involves  adding  (2.12)  to  (2.28).  This 
process  expands  the  equation  giving 

M- 1 

rnp{t)  =  aMt  -Tr-  mTr)d2l/Aj2,(n4+A)e#.  (2.29) 

m= 0 


At  this  point,  the  signal  is  at  ‘D’  in  Fig.  2.2.  The  received  signal  is  now 
ready  to  be  Matched  Filtered  (MF).  This  process  introduces  the  Time  Frequency 
Auto  Correlation  Function  (TFACF).  The  MF  process  is  the  convolution  of  the 
transmitted  signal  and  the  received  signal, 


r  nP(f)h{r 


t)dt. 


(2.30) 


where  x  represents  the  MF  output  and  h(t)  is  the  matched  filter  impulse  response. 
Since  this  model  match  filters  pulse-by-pulse,  the  impulse  response  is  simply  the  com¬ 
plex  conjugate  of  the  time  reversed  transmitted  signal,  h{t )  =  u*(—t).  Substituting 
the  received  signal  and  u*(—t)  into  (2.30)  results  in 


poo  M  1 

XnP{t)=  /  ar«(t  -  TR  -  -  r)(it.  (2.31) 

m= 0 


Pulling  the  terms  not  dependant  on  r  out  of  the  integral  gives 


M—l 


x(t)  =  arej2^n^+pd^e^  /  £  u(t  -  TR  ~  mTr)u*(t  -  r)ej2nfdtdt.  (2.32) 


m= 0 
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Adjusting  the  notation  of  the  equation,  the  summation  is  removed  from  the  equation 
and  noted  in  the  subscript. 


/»oo 

XmnP(t)  =  are^n^+p^e^  /  u(r  -  TR  -  mTr)u\r  -  t)e^Tdr  .  (2.33) 


This  underlined  section  is  the  TFACF  labeled  <f>. 

The  next  couple  of  steps  focus  on  the  TFACF.  First  to  simplify  the  TFACF, 
a  change  of  variables  is  performed  with  the  integral.  This  change  of  variables  is 
completed  by  first  letting  (3  =  t  —  TR  —  mTr  which  results  in  t  =  f3  +  TR  +  mTr  and 
=  1,  modifying  the  TFACF  expression  to 


®{T,fd)=  u(f3)u*(/3  +  TR  +  mTr-r)e 


j2nfd(0+TR  +  mTr) 


d(3. 


(2.34) 


The  bracketed  section  of  the  exponential  has  no  dependance  on  (3.  Therefore,  these 
exponentials  are  pulled  out  front  of  the  integral  and  incorporated  into  the  rest  of  the 
MF  expression.  This  leaves  the  expression,  coming  of  the  MF  at  ‘E’  in  Fig.  2.2,  for 
the  TFACF  as 


®{T,fd)=  /  u((3)u*((3  +  TR  +  7nTr-T)ej27TfdPd(3.  (2.35) 


Reviewing  the  MF  expression  (2.33)  after  simplifying  the  TFACF,  the  overall 
expression  becomes 


•v  —  n  JM^+P&z)Ji>J27rfdTR  j2irfdmTr 

Amnp\L )  u/ro  o  o  p 


u((3)u*((3  +  Tr  +  mTr  —  r)e:,27r^d/3d/3, 


(2.36) 


where  the  underlined  exponential  has  the  /^Tr  term.  As  defined  earlier,  f^Tr  equals 
the  normalized  Doppler  frequency  u>.  Therefore,  rewriting  this  equation  with  the  uj 
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notation  and  moving  some  terms  around  results  in 


,jipej2nfdTR  ej2ir(n&x+p&z) ej2nmw 


L 


■OC 


u(/3)u*(/3  +  Tr  +  mTr  -  r)ej2nfdf)df3.  (2.37) 


—  OO 


Examining  the  bracketed  portion,  these  terms  are  all  constant  terms  that  have  no 
dependance  on  time  or  pulse  number.  Therefore,  all  of  these  variables  are  rolled  into 
one  variable  at  which  provides  a  single  complex  amplitude  term  corresponding  to  a 
point  scatter.  This  simplification  leaves  the  MF  expression  as 


'OO 


L 


J2Tr(ni9x+p&z)  ej2iTmuj 


u(P)u*(/3  +  Tr  +  mTr  -  T)ej2*fd0df3.  (2.38) 


Xmnp(t)  —  CpC' 


— OO 


!  The  radar  model  will  be  continued  to  be  built  from  this  point  further 
implementing  key  assumptions.  However,  this  point  will  be  revisited  later  in 
Chapter  three  as  the  starting  point  for  exploring  the  effects  of  the  TFACF 
on  the  radar  model. 

At  this  point,  two  assumptions  are  made  within  the  TFACF.  These  assump¬ 
tions  are  that  the  received  signal  is  matched  perfect  in  range  and  in  Doppler.  The  first 
assumption  is  that  the  target  range  cell  corresponds  to  a  range  delay  r  =  TR  +  mTr. 
Breaking  this  assumption  down,  there  are  really  two  parts  to  this  assumption.  The 
first  part  is  that  there  is  a  target  present  at  this  point  in  time  and  space.  This 
assumption  will  still  be  held  through  out  this  research.  The  second  part  to  this 
assumption  is  that  when  matching  this  location,  the  range  cell  matches  perfectly. 
For  the  research,  this  assumption  is  removed  allowing  for  the  range  cell  to  vary  |  a 
chip  width.  Implementing  the  assumption  that  the  pulses  are  matched  perfectly  in 
range  into  (2.38),  the  expression  becomes 


(2.39) 
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The  second  assumption  is  that  the  transmitted  waveform  is  Doppler  tolerant. 
In  effect,  this  assumption  removes  the  Doppler  shift  exponent  from  the  TFACF 
integral,  reducing  the  integral  to 


u(/3)u*  (/3)d/3 


1. 


(2.40) 


With  the  energy  being  defined  in  (2.7)  equaling  one,  this  integral  also  evaluates  to 
one.  With  the  TFACF  equaling  one,  there  is  no  effect  on  the  MF  in  the  radar  model. 

These  assumptions  make  the  generic  mnpth  received  matched  filtered  signal 
expression  for  a  specific  range  cell  k  (target  location)  to  be 


\mnpk 


—  atej2Tv(n$x+p'&z)  ej2irmuj 


(2.41) 


Equation  (2.41)  is  the  final  expression  for  the  matched  filter  output  which 
contains  the  terms  dx,  'dz,  and  u.  These  terms  contain  information  in  space  and 
time  about  the  returned  pulse.  An  examination  of  the  role  that  these  terms  play  is 
now  necessary. 

2-4  Steering  Vector 

Equation  (2.41)  describes  the  received  signal  in  space  and  time  at  a  particular 
range  cell.  A  rearrangement  of  this  equation  into  a  linear  algebra  framework  simpli¬ 
fies  simulations.  For  insight  into  how  the  linear  algebra  will  help,  (2.41)  is  rewritten 
as 

Xk  =  atej2nl(0:N-l)H#x+(0:P-l)H#4ej2n(0:M-l)Hio  (2.42) 

where  H  represents  the  Hcrmitian  operation.  Equation  (2.42)  results  in  a  NMP  x  1 
vector  containing  all  mnp  returns  for  range  cell  k. 
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Each  vector  within  (2.42)  is  a  different  steering  vector.  These  three  steering 
vectors  describe  the  phase  progression  across  azimuth  channel  elements, a(i9x), 


a  (i9x)  = 


0j2ivdx  j2n(2)tix  .  .  .  j2ir(N-l)i}x 


elevation  channel  elements,  e(th), 


e(tfz)  = 


0j2ivdz  j2n(2)dz  .  .  .  _j27r(JV— 1)0* 


and  across  pulses,  b(u;), 


bM  = 


ei2nu  ej 2tt(2)u  .  .  .  pj2Tr(N-l)u 


(2.43) 


(2.44) 


(2.45) 


Combining  the  steering  vectors  is  done  using  the  Kronecker  product  to  form 
the  vector  in  (2.42).  This  new  vector  describing  the  phase  progression  from  element- 
to-element  and  pulse-to-pulse  is  the  steering  vector,  v, 

v  =  e(i9z)  0  b(u>)  0  a($x).  (2.46) 

Placing  the  steering  vector  into  (2.42),  the  new  space-time  snapshot  describing  each 
element  and  pulse  return  at  range  cell  k  is 

Xk  =  otte(tfz)  0  b(u))  ®  a(i9x).  (2.47) 


2.5  Element  Pattern 

To  create  the  element  pattern,  an  assumption  is  made  that  all  the  elements 
are  identical.  Identical  elements  means  each  element  has  the  same  voltage  pattern 
/(6I,0)  and  power  pattern  g(9,4 >),  where  the  power  pattern  is  equal  to  the  voltage 
pattern  squared. 
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To  create  the  element  pattern,  a  simple  cosine  pattern  in  azimuth  and  elevation 
is  used.  This  pattern  creates  a  complete  spherical  shape  surrounding  the  platform. 
Since  the  target  of  interest  is  located  in  front  of  the  radar,  there  is  no  need  to  transmit 
or  model  the  pattern  behind  the  radar.  Yet,  no  antenna  is  perfect  and  some  energy 
is  still  effectively  transmitted  behind  the  radar.  Therefore,  a  backlobe  attenuation 
factor  be  is  incorporated  in  the  element  pattern  to  represent  energy  transmitted 
behind  the  radar.  This  attenuation  factor  is  typically  modelled  30-dB  lower  than 
that  of  the  mainlobe  representing  the  concept  that  the  energy  from  the  antenna  is 
transmitted  forward  at  the  target  and  not  behind  the  radar  where  there  is  no  target. 

The  element  pattern  and  the  array  pattern  are  needed  to  be  able  to  beamform 
the  antenna.  Since  the  element  pattern  was  just  found,  the  next  step  is  to  find  the 
array  pattern. 

2.6  Array  Pattern 

The  array  pattern  uses  the  element  pattern  generated  to  form  the  antenna 
pattern  by  coherently  adding  the  signal  transmitted  from  each  element.  This  co¬ 
herent  addition  creates  the  antenna  beam.  To  steer  this  beam,  relative  time  delays 
between  elements  create  phase  delays.  These  time  delays  were  already  found  and 
are  expressed  in  the  steering  vectors  a($x)  and  e(tb)  These  phase  delays  steer  the 
beam  in  the  elevation  and  azimuth.  Upon  receive,  the  returns  are  correctly  scaled 
depending  on  how  the  antenna  beam  illuminates  them.  Therefore,  the  phase  weights 
from  the  transmit  are  applied  to  the  received  returns. 

When  creating  the  array  pattern,  only  two  of  the  steering  vectors  are  needed 
(e($~)  and  a(i?x)).  This  need  for  only  using  the  vectors  e(i9z)  and  a(i?x)  is  due  to 
the  fact  that  the  pattern  only  exists  within  the  spatial  domains.  Since  the  steering 
vectors  describe  the  phase  progression  across  each  element  in  the  array,  creating  a 
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weighted  array  factor,  W(6,<p),  for  all  elements  in  the  array  is  simply  completed  by 


ty(0,0)  =  [e(75g®  a(tf*)]Tl.  (2.48) 

Rewriting  (2.48)  in  a  different  form, 

p- 1  JV-l 

W{0,  0)  =  ej2np^ej2nn^.  (2.49) 

p= 0  n= 0 


In  an  effort  to  simplify  (2.49),  the  partial  sum  of  a  geometric  series  (X^=o  a<?  = 
)  is  applied.  Applying  the  geometric  series  is  allowed  since  there  is  a  finite 
number  of  array  elements.  This  application  results  in 


IR(M)  = 


X  —  gj27rPi?2  ^  _  ej2nN0x 


(2.50) 


Examining  (2.50),  simple  factoring  from  the  exponentials  results  in 


IE(M)  = 


e 3’nP’dz  /  g-jnPtiz  _  ejKPtiz\  ejTrNtix  /  e-jnNtfx  _  ejirNdx 


gjirtiz 


e-j-rrtfz  —  eiWz 


g  jirdx 


g  —  g  j'JT'&x 


(2.51) 


At  this  point,  the  exponentials  contained  within  the  parenthesis  are  recognized  as 
sinusoids.  Rewriting  (2.51)  with  the  sinusoids  gives 


W(Q  M  =  eMP-i)OzeMN-i)Ox  f  sin(7TP#2)  sm(rrmx) 

\  sin(7ri?2)  sin(7n?x) 


(2.52) 


Substituting  in  the  values  of  i9x  and  i9z  reveals  the  angles  that  occur  within  the  array 
pattern, 


-  (p-i)  sin(0)  cos(6»)  sin (0) 


W(6,(j))  =  e3  x°  ^  ej  A0 


sin(^P  sin  0)  sin(z^hAr  cos  8  sin  0) 


sin(pA  sin  0)  sin(2^$-  cos  9  sin  0) 


(2.53) 


ad 
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Equation  (2.53)  is  the  final  equation  for  the  array  factor.  However,  there  are 
three  special  cases  that  must  be  addressed  before  continuing.  The  ‘sin’  term  raises 
concern.  When  either  0  or  9  is  at  a  value  of  zero,  the  term  sin  equals  zero.  This 
results  in  a  [j  case. 

The  first  case  dealt  with  is  when  0  =  0.  When  0  equals  zero,  the  azimuth 
portion  of  the  array  equation  goes  to  zero.  Using  L’Hopital’s  rule,  the  derivative  of 
the  numerator  and  the  denominator  is  taken  and  the  expression  is  then  evaluated  at 
0  =  0, 


W{0, 0)  = 


(1)  (sin(^Psind) 


sin  (  ^  sin  9 


A  |0  =  cos  cos 9  cos 0J  ^ N  cos 9  cost 


ilosin  (ft  cos 0 cos 0) 


.  (2.54) 


Evaluating  the  above  expression,  the  final  response  for  the  antenna  pattern  at  0  =  0 

is 


W(0,  0)  =  NejEt(p~1)s'me 


sin  (  P  sin  9 
sin  (  ^  )  sin  9 


(2.55) 


The  second  case  dealt  with  is  when  9  =  0.  Again,  using  L’Hopital’s  rule  the 
response  for  the  antenna  pattern  is  similar  to  that  in  (2.55), 


W(O,0)  =  PejlIt{N-1)sin6 


sin  (  ^  J  N  sin  6 
sin  (  ^  )  sin  9 


(2.56) 


The  third  and  last  special  case  that  needs  to  be  addressed  is  when  both  0  and 
9  are  equal  to  zero.  Using  the  same  method  outlined  above,  the  weighting  factor 
becomes 

W (0,  0)  =  NP.  (2.57) 
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The  array  factor  is  now  modeled  accurately.  However,  it  is  symmetrical  in  the 
complete  sphere.  The  array  factor  is  multiplied  by  the  element  pattern  to  account 
for  the  back  lobe  attenuation.  The  product  of  this  union  provides  for  an  accurate 
transmit  antenna  beam  produced  by  the  radar. 

The  antenna  is  now  modeled  accurately.  This  accurate  model  is  necessary  to 
implement  later  with  the  clutter  model  allowing  for  proper  scaling  of  clutter  returns. 
However,  the  next  logical  step  is  to  examine  the  effects  that  thermal  noise  produces 
within  the  receiver  elements. 

2. 7  Thermal  Noise 

Thermal  noise  or  Additive  White  Gaussian  Noise,  is  present  at  every  array 
element.  This  noise  is  produced  naturally  within  the  elements  themselves  and  must 
be  accounted  for  within  the  model.  Calculating  thermal  noise  depends  on  the  radar 
bandwidth  and  the  standard  (operating)  temperature.  Noise  power,  typically  repre¬ 
sented  as  a2,  is  equal  to  the  product  of  the  system  bandwidth  and  the  noise  power 
spectral  density  N0.  Na  is  a  calculated  from  Boltzman’s  constant  k  and  the  standard 
temperature  Ta. 

Since  thermal  noise  is  unique  to  each  separate  element,  the  statistics  of  the 
noise  in  the  spatial  domain  are  mutually  uncorrelated.  The  temporal  domain  is  also 
mutually  uncorrelated  since  it  is  white  noise.  This  uncorrelation  results  in  a  space- 
time  snapshot  xn  describing  thermal  noise.  Finding  the  covariance  matrix  for  the 
thermal  noise  is  completed  using  xn  and  taking  the  expected  value  of  the  its  outer 
product.  However,  this  matrix  is  simplified  due  to  the  correlation  along  the  diagonal 
of  the  matrix.  Everywhere  else  in  the  matrix,  the  values  are  zero.  Therefore,  the 
thermal  noise  covariance  matrix,  Rn,  simplifies  to 

R-n  =  MNPi  (2.58) 
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where  I mnp  is  an  Identity  matrix. 

At  this  point,  the  model  is  almost  complete.  One  factor  not  incorporated  into 
the  model  still  needs  examining.  When  the  radar  transmits  a  signal,  signals  are 
returned  from  more  than  just  the  point  scatterer  target.  Returns  from  the  ground, 
sea  and  atmosphere  are  also  received  by  the  radar.  These  returns  are  called  clutter. 


2. 8  Clutter 


The  clutter  model  used  is  for  this  model  is  one  that  is  generally  accepted  within 
the  radar  community.  Sources  of  clutter  are  limited  to  those  produced  by  the  Earth’s 
surface  (land  and  sea).  The  effects  of  clutter  from  the  sea  and  ground  are  rather 
large  when  compared  to  the  clutter  from  the  atmosphere.  The  Earth  is  assumed 
perfectly  spherical  with  a  |  effective  radius  of  the  Earth. 

A  limitation  of  the  radar  model  is  only  inspecting  a  single  range  cell  at  a  given 
moment.  The  range  cell  exists  in  a  complete  circle  encompassing  the  aircraft.  This 
ring  is  called  the  clutter  ring.  Each  ring  is  divided  into  Nc  patches.  The  clutter 
model  is  a  summation  of  all  Nc  patches  creating  the  clutter  ring  for  a  specific  range 
cell.  Typically,  the  range  cell  of  interest  is  the  target  location. 


Defining  each  clutter  patch  is  the  beginning.  Each  patch  has  a  distinguished 
reflectivity  that  is  dependant  on  the  elevation  and  azimuth  angles.  To  begin  with, 
the  elevation  angle  is  calculated  9C.  Referencing  Fig.  2.3,  9C  is  found  using  the  law 
of  cosines 


9r  = 


.  -i  R2C  +  ha(ha  +  2ae) 
sm 


(2.59) 


2  Rc(ae  +  ha ) 

where  Rc  is  the  range  to  the  clutter  patch,  ha  is  the  height  of  the  aircraft,  and  ae  is 
the  effective  radius  of  the  earth. 


The  next  step  is  to  find  the  grazing  angle.  The  grazing  angle  is  defined  by  the 
tangential  line  from  the  clutter  patch  to  the  line  from  the  patch  to  the  aircraft,  as 
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Figure  2.3:  This  represents  the  clutter  geometry.  The  view  is 
from  the  side  providing  insight  into  the  calculation  of  the  grazing 
and  clutter  angles. 


seen  in  Fig.  2.3,  resulting  in  a  expression  for  the  grazing  angle  as 


,  .  _i  R2C 

-0C  =  —  Sill  — 


^a(^a  T  2ae) 

2  Rcae 


(2.60) 


Every  radar  contends  with  return  ambiguities.  This  model  incorporates  any 
possible  clutter  ring  ambiguities  that  occur  if  the  horizon  range  extends  greater  than 
the  unambiguous  range  of  the  radar  where  the  unambiguous  range  is  defined  as 
Ru  =  [5].  The  horizon  range  is  found  at  the  point  where  the  grazing  angle  equals 

zero, 

Rh  =  \/K  +  2  Haae.  (2.61) 

The  number  of  range  ambiguities  is  then  found  by 

Nr  =  (2.62) 

If  Rh  >  RcK,  Nr  is  rounded  up,  otherwise  Nr  is  rounded  down. 

At  this  point,  a  3-D  space-time  snapshot  of  the  clutter  for  a  specific  range  cell 
k  is  expressed  as 


vy-i Ay-i 

Xc  =  zL  aueu^z)  ®  ®  (2.63) 

i= 0  i=0 
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This  expression  characterizes  returns  for  every  clutter  patch  around  all  Nr  rings  for 
a  specific  range  cell  k.  However,  the  random  clutter  patch  amplitude  an  is  unknown. 

Finding  an  starts  with  determining  the  effective  RCS  of  the  clutter  patch, 
an-  The  RCS  patch  is  found  by  multiplying  the  area  of  the  patch  by  the  ground 
reflectivity,  a0(9i,(pi ),  for  the  patch, 


<?u  =  a0(9i,  ^RiA^ARsecipi. 


(2.64) 


The  ground  reflectivity  is  found  using  a  constant  gamma  model  for  clutter  [5]  and  is 
defined  as 

cro(0  (j)i)  =  7sin^i.  (2.65) 


The  area  of  the  patch  depends  on  the  range  Ri  as  well  as  the  grazing  angle  tjjr  to  the 
patch.  Acf),  defined  as  |7.  is  the  angular  extent  of  the  patch  while  A R  is  the  range 
resolution,  defined  as 

With  the  clutter  patch  RCS,  an  element  Clutter  to  Noise  Ratio  (CNR)  is 
found  using  the  radar  range  equation.  Finding  the  CNR  requires  array  and  element 
patterns.  This  need  for  the  patterns  is  why  time  was  spent  earlier  developing  those 
section  of  the  model.  The  CNR  for  a  single  element  is  defined  as 


PtGt(9i,  4>i)g(9j,  4>i)\20an 
(■ iTT)3N0BLsRf 


(2.66) 


where  Pt  is  the  power  transmitted,  Ls  is  the  loss  in  the  system,  A  is  the  wavelength, 
Gt(9ij(pi )  is  the  array  pattern,  and  g(9i,(j)i )  is  the  element  pattern. 

Now  the  random  amplitude  can  be  defined  as 


otu  =  a2^n- 


(2.67) 
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All  this  development  enables  a  3D  space-time  snapshot  to  be  rewritten  as 


JVr- 1  JVC— 1 

Xc  =  ^22  X  ®  (2.68) 

i= 0  1=0 

In  order  to  find  the  clutter  covariance  matrix,  the  expected  value  of  the  outer 
product  of  (2.68)  is  performed 


=  B{XoX?}. 


(2.69) 


which  is  rewritten  as 
Nr-  1  Nc—1 

Rc  =  22  X  °'2^«ea(^)e«  (Pz)  ®  hii(^)hu  (w)  ®  au(0x)a?[ (2.70) 

i= 0  1=0 

The  radar  model  is  now  completely  built  covering  the  radar  from  transmission 
of  the  signal  to  the  processing  of  the  received  signal.  An  acknowledgement  is  made 
that  outside  interferers  such  as  clutter  returns  and  thermal  noise  in  the  receivers  effect 
the  detection  of  targets.  Therefore,  these  interferes  are  modeled  and  incorporated 
with  in  the  radar  model.  At  this  point,  the  research  continues  onto  looking  at  the 
effects  that  the  TFACF  will  have  on  the  model. 
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III.  Time  Frequency  Auto  Correlation  Function  Developement 


The  Ambiguity  Function  (AF)  is  a  tool  to  analyze  ambiguities  that  occur  within 
waveforms.  The  AF  is  simply  the  magnitude  of  the  Time  Frequency  Auto 
Correlation  Function  (TFACF),  the  matched  filter  output  of  the  signal,  squared. 
This  matched  filter  produces  an  autocorrelation  of  the  signal.  However,  there  is  a 
Doppler  shift  that  is  incorporated  in  the  returned  signal.  This  shift  results  in  a 
cross-correlation  producing  the  overall  TFACF. 

3.1  Generic  Waveform  Model 

Looking  back  at  chapter  two,  an  assumption  was  made  to  approximate  the 
TFACF  in  (2.38)  to  one.  At  this  point,  the  expression  for  the  TFACF  is  needed  for 
placement  within  the  model. 

The  first  step  in  building  the  generic  TFACF  is  the  the  transmit  signal  given 
by 

M—l  Q-l 

»«)  =  ££  u{t  -  mTr  -  qTc)e^m\  (3.1) 

m=0  q= 0 

where  ipmq  is  the  random  starting  phase  of  the  signal,  Tr  is  the  PRI,  Tc  is  the  chip 
width  within  the  pulse,  and  Q  is  the  number  of  chips  within  each  pulse.  The  received 
signal  is  modeled  as  the  transmit  signal  with  a  time  delay  to  the  target  as  well  as  a 
Doppler  shift  on  the  signal.  This  received  signal  is  written  as 

r(t)  =  s(t  -  TR)ej27rfd(t-TR)  (3.2) 

where  TR  is  the  round  trip  time  to  the  point  scatterer. 
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Since  the  TFACF  is  the  section  under  examination,  a  closer  look  at  it  is  desired. 
The  TFACF  is  defined  as 


HrJd) 


r(t)h{r 


t)dt 


(3.3) 


where  h(t )  is  the  matched  filter  impulse  response.  Since  each  pulse  is  matched  filtered 
upon  receive,  the  matched  filter  impulse  response  is  considered  the  conjugate  of  the 
transmitted  signal,  h(t)  =  s*(—t).  Substituting  into  the  TFACF  for  r(t)  and  h(t) 
results  in 

/OO 

s(t-TR)eP2^s*(t~r)dt.  (3.4) 

-OO 

At  this  point,  Tr  and  r  are  providing  the  same  delay  action  on  the  equation.  So  the 
TFACF  equation  simplifies  to 


/OO 

s(t)ej27rfdts*(t-r  +  TR)dt.  (3.5) 

-OO 

Therefore,  inserting  the  transmit  signal  into  the  TFACF  results  in 


/»oo  M  1  Q  1 

$(t  -  Tr,  fd)=  Y.Y.  u (*  ~  mTr  ~  QTc)ej^qej2^ 

m= 0  q= 0 

M- 1  Q—l 

EE  u(t  —  aTr  —  bTc  —  t  +  TR)e  ^ abdt ,  (3.6) 

a= 0  6=0 

where  a  and  b  represent  the  returned  signal  pulse  and  chip  numbers,  respectively. 
This  is  required  to  account  for  offsets  between  the  transmit  and  received  signals.  At 
this  point,  an  assumption  is  made  that  there  are  no  multiple  time  around  returns. 
Therefore,  this  assumption  indicates  that  the  pulses  are  matched  (n  =  a),  allowing 
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one  of  the  summations  to  collapse  leaving 


M—  1  /»oo  Q  1 

Hr  -  Tr,  fd)  =  £  I  ~  mTr  ~  <W*"a 

m= 0  q= 0 

Q- 1 

u(t  -  mTr  -bTc-r  +  TR)e~j^bej2nfdtdt.  (3.7) 

6=0 

Another  assumption  is  made  at  this  point  regarding  the  TFACF.  A  limitation 
is  placed  on  the  shift  of  the  signal  as  it  is  received.  The  pulse  will  not  be  more  than 
a  half  chip  width  off,  |r  <  T-R~Tc  | .  This  assumption  allows  for  the  target  SNR  to 
remain  constant.  With  this  limitation,  there  are  three  different  cases  that  must  be 
considered.  These  cases  are  when  r  —  Tr  <  0,  r  —  TR  —  0,  and  r  —  TR  >  0.  At  this 
point,  a  deeper  examination  of  these  three  cases  is  completed. 

3.2  Examination  of  the  Time  Frequency  Auto  Correlation  Function 

There  are  three  cases  that  arise  due  to  the  limitation  placed  on  the  pulse  being 
less  than  one  half  of  a  chip  off,  r  —  TR  <  0,  r  —  TR  —  0,  and  r  —  TR  >  0. 

3.2.1  Less  than  Case  (t  —  TR  <  0).  This  case  occurs  when  the  received 
signal  is  early.  Figure  3.1  illustrates  how  the  transmitted  pulses  corresponds  to  the 
received  pulse  and  thus  how  the  chips  within  the  pulse  line  up.  Referencing  Fig  3.1 
assists  in  determining  the  limits  of  integration. 

Since  there  are  two  different  phases  from  the  received  signal  that  are  being 
matched  to  the  transmitted  signal,  two  different  integrals  are  required.  The  different 
phases  occur  due  to  the  phase  cancellation  of  the  perfectly  matched  pulse  and  the 
lack  of  phase  cancellation  with  the  offset  matched  pulse.  Referencing  Fig  3.1,  the 
limits  of  integration  are  determined. 

The  perfectly  matched  signal  is  mTr  +  qTc  for  the  lower  limit  and  the  upper 
limit  is  mTr  +  (q  +  1)TC  +  r  —  TR.  The  limits  for  the  signal  that  has  offsetting  phases 
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Figure  3.1:  The  case  is  when  the  received  chip  is  matched 

early.  The  phase  offset  are  shown  by  a  cross-hatch  while  the 
equal  phase  cases  are  blank. 


are  mTr  +  (q  +  1  )TC  +  r  —  TR  for  the  lower  limit  and  mTr  +  (q  +  1  )TC  for  the  upper 
limit.  Therefore,  the  TFACF  si  now  written  as 


M-lQ-l  „mTr  +  (g+l)Tc  +  r_Ti? 

<f>(r-T«<0,/d)  =  EE/  „  Q—j'&ma ^j27T fdt ^  | 

m=0  <?=0  JmTr+qTc 

/*mTr-\-(q-\-l)Tc 

/  eMmqe-i1>m,q+lej^fdtdt  (3.8) 

</mTT'+(q,+l)Tc+T-T^ 

The  phases  in  the  first  integral  cancel  each  other  out  leaving 


M— IQ-1  „mTr+(q+i)Tc+T_T 

EE/  _  ej2nfdtdt+ 

m= 0  q=0  J  ™Tr+qTc 


rmTr-\-(q-\-l)Tc 


'  mTr-\-(q-\-l)Tc-\-T—Tji 


e3^mq  e  3^m.iq+lei^fdtdt  (3.9) 


The  solution  to  the  first  integral  is 


aj2irfdt 


jZnfd 


mTr-\-(q-\-l-)Tc-\-T—Tji 


mTr-\-qTc 


(3.10) 


Evaluating  this  solution  by  inserting  the  limits  of  integration  in  for  t  gives 


jZnfd 


J2irfd(m.Tr+(q+l)Tc+T-TR)  _  j2irfd(mTr+qTc 


(3.11) 


Pulling  out  the  common  terms  within  the  expression  leaves 


J^fdqTc  j2-KfdmTr  \pj2irfd(Tc+T-TR) 


jZnfd 


1] 


(3.12) 


The  exponential  term  containing  Tr  is  pulled  outside  of  the  summation  over 
Q  since  it  is  dependant  on  M.  This  term  occurs  within  the  second  integral  as  well. 
Since  the  Tr  term  is  dependant  only  on  M,  it  forms  the  steering  vector  b(o;). 
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Looking  back  at  (3.12),  a  exponential  term  is  pulled  out  of  the  brackets  giving 


ej2nfdqTc 

j2vr  fd 


ej^fd(Tc+T-TR) 


\^e3^fd{Tc+r-TR)  _  e-j-rrfd(Tc+T-TR 


(3.13) 


At  this  point,  recognition  is  made  that  this  form  is  similar  to  Euler’s  identity.  The 
numerator  and  denominator  are  multiplied  by  (Tc  +  r  —  TR)  to  complete  Euler’s 
identity  allowing  the  equation  to  simplify  to 


(Tc  +  r  -  TR)ej2^dqTce^fd^+T-T^smc[fd(Tc  +  r  -  TR)}.  (3.14) 


Equation  (3.14)  is  the  final  expression  for  the  Erst  integral.  Now  the  second 
integral  needs  to  be  solved, 


r'mTr+(q-\-l)Tc 


e3^mqe-j'l/jm,q+lej2'Kfdt(lj.' 


mTr+(q+l)Tc-\-T  Tr 


(3.15) 


The  phase  exponentials  are  combined  together.  Since  they  have  no  dependance  on 
t,  the  phase  exponentials  are  pulled  outside  the  integral  resulting  in 


mq  iftm, q-\- 1 ) 


f>mTr+(q+l)Tc 


ej2nfdtdt 


>  mTr -\- ( q 1 )  Tc  + t — Tr 


(3.16) 


The  answer  for  the  solved  the  integral  is 


o  i^Pmq  '4)m,q+ 1) _ fdt 


jZnfd 


inTr  ~\~  ( q + 1 )  Tc 


mTr+(q+l)Tc-\-T—TR 


(3.17) 


Evaluating  the  solution  with  the  limits  of  integration,  the  expression  becomes 


0(3pmq  —  'lpm,q+ 1)  _ 


J'2vr  fd 


fd{mTr+{q+l)Tc)  _  ej2irfd(mTr+(q+l)Tc+T-TR) j 


(3.18) 
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Pulling  out  the  common  terms  from  the  expression  gives 


0j2nfdmTr(tl)mq—iprn,q+1) 


0j2nfd{q+l)Tc 


jZnfd 


[1 


J2-Kfd{T-TR 


>] 


(3.19) 


Again,  the  exponential  term  containing  the  mTr  is  pulled  out  ahead  of  the  summation 
over  Q  due  to  it  being  folded  into  the  steering  vector  b(u>).  Euler’s  identity  is  found 
by  pulling  a  term  out  from  the  brackets  resulting  in 


_eW  mq  "0m, g+ 1) 


ej27r/d(g+l)Tc 

j2vr  fd 


ejirfd{T-TR)  [ejTrfd(T-TR)  _  e~jirfd(T-TR)^ 


(3.20) 


Completing  the  equation  to  form  a  proper  Euler’s  identity,  the  numerator  and  de¬ 
nominator  are  multiplied  by  r  —  TR, 


-  e^mq-^m,q+ i)(r  _  TR)ej2nfdiq+1)TceMT-TR]smc[fd(r  -  TR)},  (3.21) 


which  is  the  final  form  for  the  second  integral. 

Combining  (3.14)  and  (3.21),  the  overall  expression  for  —  y  <  r  —  TR  <  0  is 

M—l  <3-1 

$(r  -  Tr  <  0,  fd)  =  ej2nfdmTr  J^(TC  +  r  -  TR)ej2nfdqTcejnfd{T°+T-Tlt) 

m= 0  q=0 

sinc[/d(Tc  +  r  -  TR)]  -  e(^-^*.«+0(T  -  TR) 

ej2nfd(q+l)TCej^fd(r-TR)s[llc^T  _  (3.22) 


3.2.2  Greater  than  Case  (r  —  TR  >  0 ).  This  case  occurs  when  the  received 
signal  appears  to  be  late.  Figure  3.2  illustrates  how  the  received  pulse  matches 
up  to  the  transmitted  pulse  as  well  as  how  the  chips  within  the  pulses  are  offset. 
Referencing  Fig  3.2  provides  insight  into  the  limits  of  integration. 

Just  as  in  the  less  than  case,  this  case  also  has  two  integrations  due  to  the  phase 
shifts.  The  first  integral  with  matching  phases,  the  lower  limit  is  mTr  +  qTc  +  r  —  TR 
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Figure  3.2:  This  case  occurs  when  the  received  chip  arrives 

late.  The  equal  phases,  which  cancel  each  other  out,  are  rep¬ 
resented  by  blank  space  while  the  offset  phases  that  interferer 
with  each  other  are  cross-hatched. 
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while  the  upper  limit  is  mTr  +  (q  +  1  )TC.  The  offset  phase  integral  has  upper  and 
lower  limits  of  mTr  +  qTc  +  r  —  Tr  and  mTr  +  qTCl  respectively,  resulting  in 


m-iq-i  mTr+(g+i)Tc 

<f>(r-T«>0,/d)  =  EE/  .  gj'hmq  e~ji’mq  gj27r/di 

m= 0  5=0  J  mTr+qTc+T-TR 

/■m,Tr+(q+ 1  )Tc+r  Tr 

/  ej^rnqe~j‘,prri'q+1ej2nfdtdt.  (3.23) 

J  mTr+(q-\-l)Tc 

By  canceling  out  the  common  phase  terms,  the  equation  becomes 


M-1Q-1  rmTr+(q+l)Tc 

EE/  .  ej2nfdtdt+ 

m= 0  <?=0  J  rnTr+qTc+T-TR 


ftmTr-\-(q-\-l)Tc-\-T—TR 


’  mTr + ( q + 1 )  Tc 


gihm^g  j^m,q+lej2ivfdt^t  (3.24) 


Examining  each  integral  separately,  the  first  integral  solution  is 


Dj27T  fd 


jZnfd 


mTr+(q+l)Tc 

mTr+qTc+T-TR 


(3.25) 


Evaluating  this  solution  with  the  limits  of  integration  results  in 


jZnfd 


aj2nfd(™Tr+(q+l)Tc)  _  j2nfd(mTr+qTc+T-TR 


(3.26) 


Pulling  out  the  common  terms  within  the  expression  leaves 


0j2-KfdmTr 


J'2tt faqTc  \pj2irfdTc  _  pj2nfd(T-TR 


jZnfd 


(3.27) 


Just  like  the  less  than  case,  the  exponential  term  containing  the  thTr  is  pulled  out 
in  front  of  the  Q  summation.  Pulling  the  r  —  Tr  term  out  of  the  brackets, 


jlnfd 


j2itfd(qTc+T-TR)  [  j2irfd(Tc-T+TR)  _ 


1], 


(3.28) 
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allows  for  mathematical  simplification  later  in  the  derivation.  In  an  effort  to  get  the 
bracketed  terms  into  a  form  that  resembles  Euler’s  identity,  the  equation  is  reworked 
giving 


ej2nfd{qTc+T-TR) 

jZnfd 


e3nfd{Tc-T+TR) 


fd.{Tc-T+TR)  _  e~jnfd(Tc-T+TR 


(3.29) 


The  numerator  and  the  denominator  by  (Tc  —  r  +  Tr),  completing  Euler’s  identity. 
This  reduction  allows  the  equation  to  reduce  to 


(Tc  —  r  +  Ti?)ej2T/d(,?Tc+r-Tii)e^/d(Tc-r+r«)sinc[/d(Tc  -  r  +  TR)\  (3.30) 


Now  that  the  solution  to  the  Erst  integral  is  obtained,  the  second  integral  needs 
to  be  solved.  The  solution  to  the  second  integral  is 


o3  (tyrnq  V’m.q+l). 


a32nfdt 


jlnfd 


mTr-\-(q+l)Tc-\-T—TR 


mTr-\-(q-\-l)Tc 


(3.31) 


Evaluating  the  solution  with  the  limits  of  integration,  the  expression  becomes 


oJ  i'llJmq  —  '4)m,q+ 1)  _ 


J'2vr/d 


^j2TT fd{mTr+(q+l)Tc+T-TR)  _  ej2n  fd{mTr+(q+l)Tc) j 


(3.32) 


Pulling  out  the  common  terms  from  the  expression  gives 


0j2-nfdmTr  j{ll>mq-‘>l>m,q+l)  . 


_  j2nfd(q+l)Tc  \  j2irfd(T—TR)  _  il 

jZnfd 


(3.33) 


The  term  containing  mTr  is  pulled  out  ahead  of  the  Q  summation  due  to  it  being 
folded  into  the  steering  vector  b(u>) 


e3  q+l)  _ gj^fd  (<J+1)TC  ^  j^fd  (t~TR)  _ 


jZnfd 


(3.34) 
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Pulling  an  exponential  out  of  the  bracketed  terms,  the  equation  becomes 


ot  (V’mq  — 1 1pm ,  q+ 1 ) 


ei27r/d(g+l)Tc 

j2  Ttfd 


e jirfd(T-TR)  ^ ej-rrfd{T-TR )  _  e~jirfd(r-TR)^ 


(3.35) 


Multiplying  the  numerator  and  denominator  by  r  —  TR  completes  the  expression  to 
match  Euler’s  identity.  Therefore,  the  expression  simplifies  to 

eMmq-i>m,q+ i)(r  _  TR)ej2nfd{g+1)TcejnfdT'smc[fd(r  -  TR)]  (3.36) 

Therefore,  combining  (3.30)  and  (3.36)  the  overall  expression  for  0  <  r  —  TR  < 

f  is 


M—l 


Q- 1 


$(t-Tr  >  0,  fd)  =  Y,  ^fdmTr  Y(T~r+TR)(r-TRy 


,j2nfd(qTc+T-TR)jirfd(Tc-T+TR) 


m= 0  g=0 

sine [fd(Tc  -  r  +  TR)]ei^-^^+l)eJ'2T/d(9+1)TceJ'7r/dr'sinc[/d(r  -  TR)]  (3.37) 


Equal  to  Case  (r  —  TR  =  0).  The  same  development  method  is 
followed  for  this  case.  However,  the  pulses  are  perfectly  matched.  Perfectly  matched 
implies  that  the  phases  of  each  chip  cancel  each  other  out  and  that  r  =  TR  as  seen 
in  Fig  3.3.  Therefore,  this  TFACF  becomes 

*{T  =  TR,fd)=  fT  ej2nfdtdt.  (3.38) 

Jo 

Since  the  energy  is  normalized  to  one,  this  equation  is  just  a  Fourier  Transform  of  a 
rect  function  over  the  interval  of  T.  Therefore,  the  expression  of  the  TFACF  for  the 
equal  to  case  is  simply  [6] 


Hr  =  TRjfd)=Tsmc(fdT).  (3.39) 
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Figure  3.3:  This  case  occurs  when  the  pulses  are  perfectly 

matched. 
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All  three  cases  (less  than,  greater  than,  and  equal  to)  are  required  to  be  evalu¬ 
ated  depending  on  the  location  of  t  —  Tr  within  the  clutter  ring,  where  \t  —  Tr\  <  y . 


3. 3  Linear  Frequency  Modulation 

Creating  the  expression  for  the  TFACF  of  the  Linear  Frequency  Modulation 
(LFM)  waveform  requires  different  mathematics.  This  requirement  is  because  LFM 
is  not  a  discrete  waveform  like  the  other  pulse  compression  waveform  techniques. 
Other  techniques  like  a  Barker  code  or  a  Gold  code  change  the  phase  from  chip  to 
chip.  LFM  is  a  frequency  ramp  and  therefore  contains  no  chips  within  the  pulse. 

Just  like  the  generic  case  with  dependance  on  chips,  LFM  has  three  different 
cases  requiring  examination  for  matching  up  pulses.  However,  due  to  the  continuous 
nature  of  the  LFM  waveform,  the  transmitted  LFM  form  is  different.  The  new 
transmitted  waveform  is 

M—l 

s{t )  =  rnTr)e^a{t-mTr)2  (3.40) 

m= 0 

where  a  is  the  chirp  rate.  Therefore,  the  TFACF  expression  becomes 


M—l  M—l 


/oo 

u(t  -  mTry^-^uit  -  yTr  -  t-Tr ) 

•OO 

e-j2na(t-mTr-r-TRrej2nfdtdt  (3.4^ 


m= 0  y= 0 


where  y  is  number  of  the  returned  pulse. 

At  this  point,  an  assumption  is  made  that  Tr  >  2 T.  This  ensures  the  trans¬ 
mitted  pulse  matches  with  only  one  other  returned  pulse.  If  this  assumption  does 
not  hold,  then  multiple  pulse  overlap  may  occur. 

3.3.1  Less  Than  Case  (r  —  Tr  <  0).  This  case  is  when  the  received  pulse 
is  early  as  seen  in  Fig  3.4.  The  LFM  case  is  more  difficult  to  determine  the  limits 
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A) 


Figure  3.4:  This  LFM  case  occurs  when  the  received  pulse  is 
early.  A  frequency  mismatch  occurs  which  effects  the  TFACF. 

of  integration.  Multiple  pulses  are  required  before  a  clear  picture  is  gained.  For  the 
less  than  case,  the  lower  limit  is  mTr  while  the  upper  limit  is  mTr  +  T  +  r  —  Tfi 

M~l  rmTr+T+T-TR 

$(r  -  Tr  <  0,  fd)  =  /  ej2na{t-mTr)*e-j2na{t-mTr-T-TR)Zej2nfdtdt 

m= 0  'mT- 

(3.42) 

The  first  step  is  to  expand  the  exponentials.  This  expansion  of  the  exponentials, 


CUlTr  ~\~T +T — Tr 


i  mTr 


0j  2'kcx.  {t2 —2tmTr  -\-mT2 ) 


e-j2-Ka{t2-2tmTr-2tT-TR+mT?+2T-TRmTr+(T-TR)2) ej2irfdt 


allows  for  simplification  of  the  integral.  Cancelation  of  like  terms  in  the  exponentials 
reduces  the  expression  to 


crnTr+T+r— Tr 


e-j2ira(-2tT-TR+2T-TRmTr+(T-TR)2)  ej2irfdtd^ 


'  mTr 


(3.44) 
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Pulling  terms  not  dependant  on  t  out  of  the  integral  gives 


pmTr+T+T—TR 

^-j2ira(2T-TRmTr+(T-TR)2)  /  e~  j  2-rr  a(-2tr -TR+ fd)t  ^ 


(3.45) 


'  mTr 


The  solution  to  the  integral  is 


?—j2irat(2T—TRmTr+(T—TR)2)  _ 


j2na(2r  -  Tr  -  fd) 


^-j2na(-2tr-TR+fd)t 


mTr+T+T—TR 


mTr 


(3.46) 


An  evaluation  of  the  solution  at  the  limits  of  integration  results  in 


:>-j2na(k2T-TRmTr+(T-TR)2)  _ 


j27ra(2r  -  TR  -  fd) 


[■ 


J2ira(2T-TR-fd)(mTr+T+T-TR)  _  j2na(2T-TR-fd)mTr 


]  (3.47) 


Pulling  the  common  terms  out  of  the  bracket  gives 


0—j2na{2r—TRmTr+{r—TR)2)  _ 


j2ira(2T  —  TR  —  fd) 

ej2Tta(2r-TR- fd)mTr  J2ira(2T-TR-fd)(T+T-TR )  _  (g  ^g) 


Pulling  out  a  exponential  from  the  bracketed  terms  gives 


:>-j2na(2T-TRm,Tr+(T-TR)2)  _ 


j2'KOi(2r  -Tr-  fd) 


0j2na(2T-TR- fd)mTr 


ljna(2T-TR-fd)(T+T-TR)  \  jrca(2T-TR-fd)(T+T-TR)  _  -j-Ka(2r-TR-fd)(T+T-Th 


]]  (3.49) 


Recognizing  the  expression  as  similar  to  Euler’s  identity,  a  multiplication  of  the 
numerator  and  denominator  is  needed.  The  equation  is  multiplied  by  T  —  r  —  TR  to 
complete  Euler’s  identity.  Thus,  providing  the  capability  to  rewrite  the  expression 
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Figure  3.5:  This  LFM  case  occurs  when  the  received  signal 
arrives  late.  This  creates  a  mismatch  with  the  frequencies. 


as 


ej2Txa{-2r-TRmTr)  e-j2-Ka{r-TR)2  ej2TxmTr[a(2r-TR-fd)} 


e^a(2  T-TR-fd)(T-r-TR){r  _  r  _  TR)sinc[a(2r  -  TR  -  fd)  (T  +  T  -  Tr)}  (3.50) 


After  simplification  of  the  equation,  the  final  expression  for  —  <  t  —  Tr  <  0  is 

M- 1 

<&(j  —  TR  <  0  fd)  —  e-32^mTrfde-j^a{r-TR)2 ejTra[2T-TR-fd)(T-r-TR) 

m= 0 

(T  -  t  -  TR)sinc[a(2r  -TR-  fd)(T  +  r  -  TR)]  (3.51) 


3.3.2  Greater  Than  Case  (r  —  TR  >  0j.  This  case  occurs  when  the  received 
pulse  is  late  as  seen  in  Fig  3.5.  For  the  greater  than  case,  the  lower  limit  is  mTr  + 
t  —  Tr  while  the  upper  limit  is  mTr  +  T, 

rmTr+T 

<F(t  -  Tr  >  0,  fd)  =  /  ej2™(t-mTr)2e-j2na(t-mTr-T-TR)2ej2„fdtdt  (3.52) 

m= 0  3  mTr+T-TR 
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Before  solving  the  integral,  an  expansion  of  the  exponentials, 


r»mTr+T 


'  mTr+T—TR 


ej2TTa(t2  -2t.mTr+mT?)  e-j2ira(t2 -2tmTr-2tT-TR+mT? +2T-TRmTr+(T-TR)2)  ej2n  fdt 

(3.53) 


is  completed  to  simplify  the  integral.  Cancelation  of  like  terms  in  the  exponential 
reduces  the  expression  to 


r»mTr+T 


'  mTr+r—Tj 


e-j2na(-2tT-TR-\-2r-TRmTr-\-(r-TR)2)ej27Tfdt^ 


R 


(3.54) 


Pulling  the  terms  not  dependant  on  t  out  from  the  integral  gives 


rmTr+T 


^-j27ra(2T-TRmTr-\-{r-TR)2)  /  e~j2ncy{-2tT-TR-\-fd)t ^  ^  55^ 


'  mTr-\-r—T} 


R 


The  solution  to  the  integral  is 


3-j27ra(2r-TflmTr+(r-rH)2)  _ 


j27ra(2r  -TR-  fd) 


j2irat(—2tT—TR+ f d)t 


mTr+T-TR 


.  (3.56) 


mTr+T 


Evaluating  the  solution  at  the  limits  of  integration  gives 


3-j27ra(2r-TRmTr+(r-Tfl)2)  _ 


j27ra(2r  -  TR  -  fd) 

[ej2ira(2T-TR-fd)(mTr+T-TR)  _  ei2ira(2r-TR- fd){mTr+T)^  (g  gy) 


Pulling  the  common  terms  out  from  the  bracketed  terms  results  in 


:1-j2na(2T-TRm,Tr+(T-TR)2)  _ 


j2ira(2T  -TR-  fd) 

ej2na(2T-TR-fd)(mTr+T-TR)  J2na(2T-TR-fd)(T-T-TR )  _  (g  gg) 
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Rewriting  the  equation  after  pulling  out  an  exponential  from  the  bracket  terms  gives 


e-j2ira(2T-TRmTr+(T-TR)2) 


ej2ira(2T-TR-fd)(mTr+T-TR) 

j2ira(2r  -TR-  fd) 


ejira(2T-TR-fd)(T-T-TR)  [ ej-Ka[2r-TR- fd)(T-r-TR )  _  e~jna(2T-TR- fd)(T-r-TR) j  (3  gg) 

Recognizing  a  form  similar  to  Euler’s  identity,  the  numerator  and  denominator  are 
multiplied  by  T  —  t  —  TR, 

e—j2Tra(2r-TRmTr+(T-TR)2)ej2na(2T-TR-fd.)(mTr+T-TR) 

eW2r-TWd(T-r-T«)(r  _  r  _  TR)sinc[a(2T  —  TR  —  f d){T  —  t  —  Tr )],  (3.60) 
completing  Euler’s  identity.  Simplifying  this  equation  to  the  final  form, 

e-j2ira(T-TR)2  ej2ir(2aT-TR+fd)T-TR 

eM2aT-TR+fd){T-T-TR)(T  _  T  _  TR)sinC[a'(2r  -  TR  -  fd)(T  -  r  -  TR)]  (3.61) 


M- 1 

Y  e-j2vmTrfde-j2wa{T-TR)2ej7cai2T-TR-fd)(-T-T-TR\T-r-TR)smc[a(2r-TR-fd)(T+r-TR)} 

m= 0 


(3.62) 


3.3.3  Equal  to  Case  (r  —  TR  =  0).  For  the  equal  case,  the  pulses  are 
perfectly  matched.  Therefore,  r  —  TR  =  0  which  makes  the  equation  reduces  to 

*{r  =  TR,fd)=  [T  C2^dtdt  (3.63) 

Jo 

Since  the  energy  in  the  transmitted  LFM  is  normalized  to  one,  this  equation  is  seen  as 
a  Fourier  Transform  of  a  rect  function  across  duration  T.  Therefore,  the  expression 
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for  the  TFACF  when  r  =  TR  is  [6] 


Hr  =  TRJd)  =  Tsmc(fdT)  (3.64) 

At  this  point,  all  the  derivations  for  this  research  have  been  completed.  The 
next  step  is  to  simulate  the  derived  TFACF’s  in  Ad  AT  LAB®  and  complete  analysis. 
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IV.  Results 


Results  from  examining  radar  model  effects  due  to  removing  the  Doppler  and 
range  mismatch  assumptions  are  provided  this  chapter.  Incorporating  a  non- 
approximated  value  of  the  Time  Frequency  Auto  Correlation  Function  (TFACF) 
into  the  model  induces  adverse  effects  on  performance.  An  analysis  is  completed  to 
evaluate  the  effects  that  the  TFACF  has  on  the  model.  This  analysis  is  completed 
by  comparing  the  conventional  model  to  the  results  from  the  scaled  model.  The 
simple  rectangular  pulse  was  used  to  ensure  that  code  used  to  model  the  TFACF 
was  correctly  completed.  After  the  results  are  confirmed,  several  other  waveforms 
are  implemented,  including  13-length  Barker,  26-length  concatenated  Barker,  two 
different  Gold  codes,  Linear  Frequency  Modulated  code,  and  two  different  sized 
Frank  codes.  Before  the  analysis  begins,  an  understanding  of  the  analysis  tools  is 
beneficial  for  drawing  conclusions. 


4.1  Ajialysis  Tools 

There  are  four  different  methods  used  to  gain  insight  into  the  differences  be¬ 
tween  the  approximated  model  and  the  scaled  model.  These  tools  help  to  determine 
any  performance  variations  that  the  radar  model  may  incur  by  removing  the  TFACF 
approximation  from  the  model.  These  techniques  operate  on  the  noise  and  interfer¬ 
ence  covariance  matrix  built  from  the  received  signal,  as  outlined  in  Chapter  2. 


4-1.1  Signal  Match.  The  Signal  Match  (SM)  technique,  also  known  as 
the  Fourier  Spectra,  determines  a  maximum  whenever  the  steering  vector  matches 
a  signal  vector  within  the  covariance  matrix  [4].  The  SM  technique  determines  the 
maximum  by  using  the  normalized  power  output,  defined  as 


p  = 

1  sm 


v^R  v 
vA  ’ 


(4.1) 
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where  v  is  the  steering  vector,  H  represents  the  hermitian  operation,  and  R  is  the 
covariance  matrix.  This  technique  is  a  two-dimensional  Fourier  Transform,  a  uniform 
linear  array  and  a  radar  having  constant  PRF  [4],  The  SM  is  a  low  resolution 
technique  providing  information  on  detection  ability  by  showing  the  strength  of  the 
clutter  returns. 

Caution  must  be  used  when  using  this  technique.  A  lack  of  resolution  may  in¬ 
troduce  biased  estimates  when  multiple  signals  are  present  which  may  cancel  weaker 
signals  [4].  A  byproduct  of  the  Fourier  transform,  the  SM  technique  produces  side- 
lobes  which  may  hide  the  true  clutter  response. 

This  low  resolution  technique  provides  an  accurate  picture  of  the  clutter  spec¬ 
trum,  yet  it  is  not  reliable  for  accurate  estimation  of  range  resolution.  To  achieve  an 
accurate  resolution  plot,  a  high  resolution  technique  is  needed,  which  provides  the 
ability  to  distinguishing  multiple  targets. 

4-1.2  Minimum  Variance  Estimator.  The  Minimum  Variance  (MV)  esti¬ 
mator  is  considered  a  high  resolution  technique.  MV  output  power  is  defined  as 

Pmv=(vh  R-Vt1.  (4.2) 

A  primary  difference  between  the  MV  and  the  SM  is  that  the  MV  does  not  provide 
an  accurate  magnitude  response,  implying  that  the  strength  of  the  clutter  returns  is 
not  accurate  on  these  plots.  What  MV  plots  provide  is  the  capability  to  distinguish 
multiple  targets  [4] .  Higher  range  resolution  allows  for  separation  of  targets  that  are 
located  close  together. 

A  substitution  for  the  MV  is  the  Maximum  Entropy  Method  (MEM),  which  is 
also  considered  a  high  resolution  technique.  However,  according  to  Klemm  [4],  the 
results  from  using  the  MEM  are  “quite  ‘spiky’  and  do  not  reflect  the  true  spectrum 
properly.  As  will  be  shown... the  MV  estimator  is  more  appropriate.” 
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A  combination  of  the  low  and  high  resolution  techniques  provide  good  insight 
into  the  characteristics  of  the  clutter  returns.  However,  a  tool  is  needed  to  understand 
the  difficulty  of  the  clutter  environment.  This  tool  is  the  Eigenvalue  decomposition. 

4-1.3  Eigenvalue  Decomposition.  To  determine  how  difficult  the  clutter 
environment  is,  an  eigenvalue  decomposition  is  performed.  The  eigenvalue  decom¬ 
position  simply  isolates  the  eigenvalues  of  the  covariance  matrix.  These  values  are 
then  sorted  in  ascending  order  and  plotted.  Looking  at  a  typical  plot,  the  greater 
the  number  of  eigenvalues  with  a  value  exceeding  that  of  the  noise  floor,  the  tougher 
the  clutter  environment  is.  A  tough  clutter  environment  increases  the  chance  that 
weaker  target  returns  are  masked  by  higher  clutter  returns. 

The  complete  picture  of  how  the  TFACF  effects  the  radar  is  almost  formed. 
However,  there  is  one  more  tool  that  is  needed  to  help  provide  the  complete  picture. 

4-1.4  Optimum  Matched  Filter  Output  SINR.  The  optimum  matched  fil¬ 
tered  output  Signal  to  Interference  plus  Noise  Ratio  (SINR)  is  a  function  of  the  noise, 
interference  covariance  matrix  and  the  steering  vector,  v  [2],  The  optimum  SINR  is 
computed  using 

SINRopt  =  a%vH  R^v  (4.3) 

where  is  the  per  pulse  Signal  to  Noise  Ratio  (SNR)  and  a2  is  the  noise  power. 
This  analysis  shows  the  performance  degradation  that  occurs  within  the  SINR. 

A  side  step  is  require  to  provide  more  insight  into  the  SINR  plots.  The  thermal 
noise  power  limit  is  determined  using  the  best  case  power  level  if  there  was  no  clutter 
returned.  The  thermal  noise  power  limit  is  given  by 

ffiimit  =  (4.4) 

which  simplifies  to 

ffilnit  =  ftV^v.  (4.5) 
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The  model  in  this  work  assumes  that  =  1.  Therefore,  (4.5)  produces  a  vector  of 
size  MNP  or  16  x  16  x  2  =  64.  Converting  the  values  of  this  vector  to  decibels  gives 

101og10(MAhP)  =  101og10(16  x  16  x  2)  =  27.09  dB.  (4.6) 

This  states  that  with  the  physical  radar  configuration  the  system  cannot  achieve  a 
better  SINR  than  27.09  dB,  the  thermal  noise  limit. 

Therefore,  knowing  the  thermal  noise  limit,  a  clearer  understanding  of  the 
effects  from  the  clutter  environment  is  gained.  Any  variation  below  the  thermal 
noise  limit  is  directly  related  to  clutter. 

The  techniques  outlined  above  provide  a  complete  picture  of  the  clutter  en¬ 
vironment.  This  picture  provides  insight  into  effects  from  removing  the  TFACF 
approximation  in  the  model. 

4-2  Simple  Rectangular  Pulse 

The  baseline  waveform  is  the  rectangular  sinusoidal  pulse  without  any  phase 
coding.  This  waveform  will  have  no  adverse  effects  due  to  Doppler  or  range  shifting. 
Therefore,  before  any  modeling  with  different  waveforms  is  completed,  the  TFACF 
scaled  model  is  compared  against  the  current  approximated  model,  ensuring  an  ac¬ 
curate  model. 

Analysis  starts  by  looking  at  TFACF  results.  The  TFACF  model  is  based 
on  the  expressions  derived  in  Chapter  3.  Plotting  the  TFACF,  a  visual  check  is 
performed,  confirming  that  the  results  from  the  derived  TFACF  model  in  Chapter  3 
matches  the  proven  TFACF  plots.  Examining  Fig  4.1,  the  plot  shows  exactly  what  is 
expected,  a  rectangular  correlated  signal  producing  a  triangle  shape  combined  with 
the  sine  function  as  seen  in  [5].  Therefore,  confidence  is  gained  with  the  derivation 
of  the  TFACF  and  that  it  is  correctly  modeled. 
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Figure  4.1:  The  Time  Frequency  Auto  Correlation  Function 
Plot  of  a  simple  rectangular  pulse. 

The  first  comparison  made  is  that  of  the  signal  match  technique.  Plotted  SM 
results  for  the  simple  rectangular  case  with  no  phase  coding  is  seen  in  Fig  4.2  (a). 
As  expected,  the  two  signals  from  the  approximated  radar  model  and  the  TFACF 
scaled  model  are  identical  to  each.  This  is  due  to  no  pulse  compression  techniques 
added  to  the  signal  which  means  there  is  no  special  phase  coding  structure  with  in 
the  pulse.  Therefore,  any  Doppler  shift  that  has  occurred  will  have  no  adverse  effects 
on  the  TFACF  model. 

Examining  the  minimum  variance  estimator  for  the  rectangular  pulse,  Fig  4.2  (b) 
shows  the  two  different  models  plotted  on  top  of  each  other.  Again,  this  plot  is  ex¬ 
pected.  Since  there  is  no  phase  coding  within  the  pulse,  the  TFACF  equals  one.  This 
makes  the  models  identical.  Hence  the  resolution  curves  are  equal  to  each  other. 

The  eigenvalue  decomposition,  Fig  4.2  (c),  shows  no  difference  between  the  two 
models.  This  is  expected  due  to  the  lack  of  pulse  compression  technique.  As  stated 
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(a) 


(b) 


(c)  (d) 


Figure  4.2:  Unmodulated  Rectangular  Pulse 

(a)  Signal  match  describing  the  clutter  spectrum. 

(b)  MVE  describing  the  resolution  capability  utilizing  this  waveform. 

(c)  Eigenvalue  decomposition  illustrating  the  rank  of  the  clutter  spectrum  and  the 
overall  difficulty  of  the  problem. 

(d)  Output  SINR  of  fully  adaptive  Space  Time  Adaptive  Processing  (STAP)  illus¬ 
trating  target  detection  capability  at  all  normalized  Doppler  frequencies. 
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above,  each  model  is  identical.  Therefore,  the  eigenvalues  in  each  covariance  matrix 
of  the  model  will  have  the  same  values. 

The  output  SINR  for  the  fully  adaptive  STAP  technique,  Fig  4.2  (d),  shows 
no  variance  between  the  two  models,  as  expected.  However,  the  effects  of  clutter 
are  shown  now.  Main  beam  clutter,  located  about  0  Hz,  degrades  the  SINR  level  by 
nearly  20  dB. 

With  the  baseline,  simple  rectangular  pulse  tested,  it  is  shown  that  the  scaled 
TFACF  model  is  derived  and  simulated  correctly.  By  duplicating  and  matching 
results  achieved  from  the  approximated  model,  confidence  is  gained  that  the  scaled 
model  is  correct.  Therefore,  additional  waveforms  are  next  examined. 

4-3  Barker  Code 

The  Barker  code  is  a  phase  coding  scheme  producing  equally  low  sidelobe  lev¬ 
els  [5].  Using  this  code  allows  for  better  resolution  along  the  range  axis.  The  problem 
with  the  Barker  code  is  that  the  maximum  code  length  is  limited  mathematically  to 
thirteen  sub-phases.  Typically,  for  pulse  compression,  a  large  code  length  is  desired 
to  increase  the  range  resolution.  The  Barker  coding  scheme  is  examined  as  a  starting 
reference  point  for  basic  pulse  compression  coding. 

4-3.1  13-Length  Barker  Code.  Before  examining  the  performance  charac¬ 
terization  metrics  obtained  by  implementing  the  Barker  code  into  the  model,  the 
TFACF  of  a  13-length  Barker  code  is  plotted  as  a  check.  This  check  ensures  that 
the  TFACF  correctly  models  the  Barker  code.  The  Barker  code  is  1111100110101  [5] 
and  is  coded  as  phase  shifts  within  each  chip  of  a  transmitted  pulse.  The  value 
of  1  is  represented  as  e ■J7r  and  0  is  expressed  as  ej0  or  0.  Figure  4.3  confirms,  by 
matching  expected  results,  that  the  TFACF  model  built  still  holds  true  for  the  13- 
length  Barker  code  sequence.  With  more  confidence  gained  in  the  TFACF  model, 
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Figure  4.3:  The  Time  Frequency  Auto  Correlation  Function 
plot  of  the  13-Length  Barker  Coded  pulse. 

the  comparison  tools  are  implemented  to  gain  insight  to  any  performance  trade-offs 
incurred. 

Starting  the  comparison  of  the  two  models  begins  with  the  first  technique,  the 
SM.  The  results  from  the  signal  match  are  shown  in  Fig  4.4  (a).  Immediately,  a 
reduction  of  3-dB  of  clutter  power  return  is  seen  which  indicates  there  is  less  clutter 
power  entering  the  system. 

The  next  step  is  to  determine  the  resolution  capability  with  the  Barker  code. 
Resolution  capabilities  is  shown  using  the  minimum  variance  estimator  technique. 
The  results  from  completing  the  MVE  is  shown  in  Fig  4.4  (b).  The  TFACF  scaled 
factor  return  is  not  as  sharp  when  compared  to  the  approximated  TFACF.  This  gain 
in  width  indicates  that  the  resolution  capability  is  degraded  due  to  adding  the  scale 
factor. 

The  difficulty  of  the  clutter  problem  is  shown  by  eigenvalue  decomposition, 
seen  in  Fig  4.4  (c).  According  to  the  plot,  the  scaled  factor  is  affected  more  by 
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(a)  (b) 


(c)  (d) 


Figure  4.4:  13-Length  Barker  Coded  Pulse 

(a)  Signal  match  describing  the  clutter  spectrum. 

(b)  MVE  describing  the  resolution  capability  utilizing  this  waveform. 

(c)  Eigenvalue  decomposition  illustrating  the  rank  of  the  clutter  spectrum  and  the 
overall  difficulty  of  the  problem. 

(d)  Output  SINR  of  fully  adaptive  Space  Time  Adaptive  Processing  (STAP)  illus¬ 
trating  target  detection  capability  at  all  normalized  Doppler  frequencies. 
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thermal  noise.  The  approximated  TFACF  noise  floor  is  —143.5  dB  while  the  scaled 
factor  TFACF  is  around  131.5  dB.  This  difference  in  thermal  noise  is  due  to  the 
change  in  bandwidth.  The  uncompressed  bandwidth  of  the  radar  is  defined  as 

B=f  =  1.25  MHZ.  (4.7) 

Pulse  compression  effectively  increases  the  bandwidth  by  dividing  T  into  equal  dura¬ 
tion  subpulses  by  the  number  or  the  code  length.  In  this  case,  the  code  length  is  13 
long.  Therefore,  the  compressed  pulse  width  is  This  makes  the  new  compressed 
bandwidth  to  be 

13 

Bc=-  =  16.25  MHZ  (4.8) 

which  is  15  MHz  larger  than  the  uncompressed  pulse  width.  Therefore,  a  high 
thermal  noise  level.  This  increase  in  thermal  noise  dominates  more  of  the  clutter 
returns  creating  an  easier  clutter  problem. 

Examining  Fig  4.4  (d),  the  optimum  SINR  values  are  shown.  As  mentioned 
above,  main  beam  clutter  drastically  effects  the  output  SINR  of  fully  adaptive  STAP. 
However,  by  increasing  the  bandwidth,  6  dB  SINR  improvement  is  gained. 

4-3.2  26-Length  Concatenated  Barker  Code.  In  an  effort  to  harness  the 
benefits  of  the  Barker  code,  yet  create  a  longer  pulse  than  the  original  13-length 
Barker  code,  a  concatenated  26-length  Barker  code  was  created.  Creating  this  con¬ 
catenated  code  was  accomplished  by  taking  the  original  13-length  code  and  attaching 
a  flipped  version  of  the  original  code  which  results  in  11111001101011010110011111. 
Again,  the  1  is  coded  as  eJ7r  and  0  is  coded  as  0.  The  TFACF  of  this  concatenated 
code  is  seen  in  Fig  4.5. 

Changing  from  the  13-length  Barker  the  26-length  concatenated  code  effects 
the  returned  clutter  power  as  seen  in  Fig  4.6  (a).  The  clutter  power  is  reduced 
7  dB  and  a  difference  of  2  dB  is  seen  between  the  two  different  Barker  codes.  This 
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Figure  4.5:  The  Time  Frequency  Auto  Correlation  Function 
Plot  of  the  26-Length  Concatenated  Barker  Coded  pulse. 

difference  shows  that  as  the  length  of  the  code  increases,  the  signal  strength  from 
the  clutter  return  slowly  degrades. 

The  resolution  capability  of  the  concatenated  Barker  code  shown  in  Fig.  4.6  (b) 
remains  wider  than  that  of  the  approximated  model.  Comparing  the  26-length  to 
that  of  the  13-lcngth,  a  change  in  resolution  is  negligible.  The  only  difference  in  the 
plots  is  that  the  floor  of  the  scale  factor  TFACF  increases  by  1  dB. 

The  new  bandwidth  from  the  26-length  concatenated  Barker  code  is 

Bc=^=  32-5  MHZ.  (4.9) 

The  compressed  pulse  bandwidth  has  increased  by  almost  30  MHZ.  The  eigenvalue 
decomposition  shows  in  Fig.  4.6  that  with  the  larger  bandwidth  of  the  concatenated 
pulse,  the  noise  floor  has  risen  to  129  dB,  effectively  making  the  clutter  environment 
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(a)  (b) 


(c)  (d) 


Figure  4.6:  26-Length  Concatenated  Barker  Coded  Pulse 

(a)  Signal  match  describing  the  clutter  spectrum. 

(b)  MVE  describing  the  resolution  capability  utilizing  this  waveform. 

(c)  Eigenvalue  decomposition  illustrating  the  rank  of  the  clutter  spectrum  and  the 
overall  difficulty  of  the  problem. 

(d)  Output  SINR  of  fully  adaptive  Space  Time  Adaptive  Processing  (STAP)  illus¬ 
trating  target  detection  capability  at  all  normalized  Doppler  frequencies. 
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less  difficult.  Comparing  the  two  Barker  code  TFACF,  the  concatenated  Barker’s 
noise  floor  increased  3  dB. 

The  optimized  SINR  output  plot,  Fig.  4.6  (d),  shows  a  8  dB  gain.  This  gain 
is  due  to  less  clutter  power  returned.  Comparing  the  difference  between  the  two 
Barker  codes,  the  concatenated  code  has  a  2  dB  greater  SINR. 


4-4  Frank  Code 

The  Frank  code  is  considered  a  polyphase  code.  A  polyphase  code  occurs  when 
the  pulse  compression  need  is  not  restricted  to  the  binary  coding  level  of  0  and  n  [5] . 
A  Frank  code  is  a  M  x  M  matrix  where  the  first  column  is  zeros  while  the  rest  of 
the  values  are  found  by  taking  the  number  of  each  row  of  the  matrix  (counting  up 
from  zero  to  M  —  1)  and  multiplying  by  a  phase  of  ||  radians  or  ^  degrees.  A 
reorganization  of  this  matrix  for  transmission  results  in  a  1  x  M 2  vector.  Therefore, 
a  single  uncompressed  pulse  transmission  looks  like  a  transmission  of  M 2  pulses. 


4-4-1  5-Length  Frank  Code.  The  first  Frank  code  examined  is  the  5-length 
Frank  code.  As  described  above,  the  length  of  the  actual  transmit  pulse  is  52  or  25 
subpulses.  The  phases  are  incremented  in 
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Therefore,  the  matrix  for  the  Frank  code  for  M  =  5  is 
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where  the  numbers  are  in  units  of  radians  or 
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where  the  numbers  have  units  of  degrees. 

Since  this  is  a  polyphase  code,  the  sidelobe  levels  are  lower  than  a  typical 
binary  code  and  the  code  tends  to  be  more  Doppler  tolerant.  The  TFACF  for  this 
code  is  shown  in  Fig  4.7.  This  plot  depicts  a  typical  Frank  code  showing  lower 
sidelobe  levels  in  the  outlying  regions. 

Beginning  the  analysis,  the  signal  match  comparison  in  Fig  4.8  (a)  is  examined. 
A  6  dB  degradation  in  clutter  power  is  achieved  with  this  Frank  code.  Comparing 
this  plot  to  that  of  the  concatenated  Barker  plots,  similarities  are  seen.  The  mag¬ 
nitude  and  overall  spectrum  appear  as  almost  identical  plots.  Fig  4.8  (b)  shows  the 
comparison  using  the  MVE  technique.  Again,  it  is  seen  that  the  approximated  reso¬ 
lution  is  better  than  the  scaled  model.  Comparison  of  this  Frank  code  to  that  of  the 
concatenated  Barker  show  negligible  differences.  A  picture  is  starting  to  form  that 
the  major  effect  on  the  model  is  that  of  increased  bandwidth,  allowing  for  greater 
influence  from  thermal  noise. 
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Figure  4.7:  The  Time  Frequency  Auto  Correlation  Function 
Plot  of  the  5-Length  Frank  Coded  pulse. 

The  compressed  bandwidth  of  this  Frank  code  is 

Bc  =  ^  =  31.25  MHZ.  (4.14) 

This  is  a  difference  of  30  MHZ  from  the  uncompressed  bandwidth.  This  difference 
is  seen  in  the  eigenvalue  decomposition  Fig  4.8  (c).  This  plot  shows  a  noise  floor  of 
129  dB  which  is  equal  to  that  of  the  concatenated  Barker  pulse.  This  correlation  is 
starting  to  show  that  the  larger  the  bandwidth,  the  easier  clutter  environment  the 
radar  deals  with  since  the  noise  floor  dominates  many  of  the  clutter  returns. 

Examining  the  output  SINR  of  fully  adaptive  STAP,  Fig  4.8  (d),  a  7  dB  gain  is 
achieved  by  implementing  the  Frank  code.  This  improvement  is  due  to  the  smaller 
clutter  power  returns  as  seen  in  the  SM  plots.  Again,  a  comparison  of  the  Frank  code 
results  to  that  of  the  concatenated  Barker  are  virtually  the  same  plot  suggesting  that 
the  biggest  factor  effecting  the  model  is  an  increase  in  thermal  noise  power. 
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(a)  (b) 


Optimized  SINR  (5-Length  Frank  Code  Pulse) 


(c)  (d) 


Figure  4.8:  5-Length  Frank  Coded  Pulse 

(a)  Signal  match  describing  the  clutter  spectrum. 

(b)  MVE  describing  the  resolution  capability  utilizing  this  waveform. 

(c)  Eigenvalue  decomposition  illustrating  the  rank  of  the  clutter  spectrum  and  the 
overall  difficulty  of  the  problem. 

(d)  Output  SINR  of  fully  adaptive  Space  Time  Adaptive  Processing  (STAP)  illus¬ 
trating  target  detection  capability  at  all  normalized  Doppler  frequencies. 
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4-4-2  10-Length  Frank  Code.  The  second  Frank  code  examined  is  the  10- 
length  Frank  code.  As  described  above,  the  length  of  the  actual  transmit  pulse  is 
102  or  100  subpulses.  The  phases  are  incremented  in 
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where  the  numbers  are  in  units  of  radians  or 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

36 

72 

108 

144 

180 

216 

252 

288 

324 

0 

72 

144 

216 

288 

0 

72 

144 

216 

288 

0 

108 

216 

324 

72 

180 

288 

36 

144 

252 

0 

144 

288 

72 

216 

0 

144 

288 

72 

216 

0 

180 

0 

180 

0 

180 

0 

180 

0 

180 

0 

216 

72 

288 

144 

0 

216 

72 

288 

144 

0 

252 

144 

36 

288 

180 

72 

324 

216 

108 

0 

288 

216 

144 

72 

0 

288 

216 

144 

72 

0 

324 

288 

252 

216 

180 

144 

108 

72 

36 

where  the  numbers  have  units  of  degrees. 

Increasing  the  Frank  code  to  a  length  10  increases  the  length  of  the  transmit 
codes  to  a  length  of  100.  This  increases  the  bandwidth  while  suppression  of  the 
sidelobes  is  greater  than  that  of  the  5-length  Frank,  as  seen  in  Fig  4.9. 

Figure  4.10  (a)  shows  the  signal  match  metric  for  the  10-length  Frank  code. 
Implementing  this  compression  scheme,  return  clutter  power  is  reduced  by  11  dB. 
This  is  the  greatest  hit  seen  thus  far.  However,  this  stands  to  reason  due  to  the 
greatest  increase  in  bandwidth. 

The  MVE  plot,  Fig  4.10  (b),  shows  that  resolution  is  lost  using  this  Frank  code. 
The  floor  has  risen  to  —10  dB  and  the  beam  has  fattened,  decreasing  the  resolution 
capability.  When  comparing  against  Fig  4.8  (b),  the  noise  floor  differs  by  2  dB. 

The  compressed  bandwidth  for  this  case  of  the  Frank  code  is 


Bc  =  ^  =  125  MHZ. 


(4.19) 
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Figure  4.9:  The  Time  Frequency  Auto  Correlation  Function 
Plot  of  the  10-Length  Frank  Coded  pulse. 

This  is  increase  in  bandwidth  by  over  120MHz.  Increasing  the  bandwidth  greatly 
simplifies  the  clutter  difficulty.  This  simplification  is  due  to  an  increase  the  thermal 
noise  level  as  seen  in  Fig  4.10  (c).  This  higher  noise  floor  results  in  fewer  clutter 
values  being  seen  in  the  model. 

By  implementing  the  10-length  Frank  code,  the  output  SINR  of  a  fully  adaptive 
STAP  process  increased  by  10  dB  as  seen  in  Fig.  4.10  (d).  The  SINR  for  the  scaled 
factor  TFACF  is  only  10  dB  lower  than  that  of  the  thermal  noise  limit,  the  level 
where  the  best  SINR  is  achieved. 

4-5  Gold  Code 

The  next  compression  scheme  examined  is  the  Gold  code.  The  Gold  code  has 
been  shown  to  provide  improvement  with  clutter  suppression,  unambiguous  range 
properties  and  resolution  [3].  Therefore,  a  31-length  Gold  Code  was  implemented  to 
see  the  effects,  if  any,  on  the  TFACF  scaled  model. 
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(a)  (b) 


(c)  (d) 


Figure  4.10:  10-Length  Frank  Coded  Pulse 

(a)  Signal  match  describing  the  clutter  spectrum. 

(b)  MVE  describing  the  resolution  capability  utilizing  this  waveform. 

(c)  Eigenvalue  decomposition  illustrating  the  rank  of  the  clutter  spectrum  and  the 
overall  difficulty  of  the  problem. 

(d)  Output  SINR  of  fully  adaptive  Space  Time  Adaptive  Processing  (STAP)  illus¬ 
trating  target  detection  capability  at  all  normalized  Doppler  frequencies. 
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Figure  4.11:  The  Time  Frequency  Auto  Correlation  Function 
Plot  of  the  First  Gold  Coded  pulse. 

4-5.1  First  Gold  Code.  The  first  Gold  code  implemented  is  [3] 

0001101111011010001111110100000.  (4.20) 

For  coding  purposes,  1  is  coded  as  a  phase  of  e ^  and  0  is  coded  as  0.  Figure  4.11 
shows  the  TFACF  of  the  Gold  code. 

Referencing  Fig  4.12  (a),  the  clutter  spectrum  of  the  gold  code  is  plotted. 
Examining  this  plot,  a  7  dB  degradation  of  the  clutter  return  power  is  gained  by 
implementing  this  compression  scheme. 

The  minimum  variance  plot  of  the  gold  code,  Fig.  4.12  (b),  indicates  that 
resolution  capabilities  are  lost  by  using  a  compressed  pulse  width. 

Looking  at  the  bandwidth  of  the  Gold  code, 

31 

Bc  =  -  =  38.75  MHZ,  (4.21) 
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(a)  (b) 


(c)  (d) 


Figure  4.12:  First  Gold  Coded  Pulse 

(a)  Signal  match  describing  the  clutter  spectrum. 

(b)  MVE  describing  the  resolution  capability  utilizing  this  waveform. 

(c)  Eigenvalue  decomposition  illustrating  the  rank  of  the  clutter  spectrum  and  the 
overall  difficulty  of  the  problem. 

(d)  Output  SINR  of  fully  adaptive  Space  Time  Adaptive  Processing  (STAP)  illus¬ 
trating  target  detection  capability  at  all  normalized  Doppler  frequencies. 
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an  increase  from  the  uncompress  pulse  bandwidth  by  more  than  37  MHZ  is  found. 
Therefore,  more  thermal  noise  is  introduced  into  the  system.  The  eigenvalue  decom¬ 
position  shows  this  increase  in  thermal  noise,  Fig.  4.12  (c).  The  noise  floor  is  at 
—  127  dB  which  suppresses  returns  from  the  clutter  environment,  making  the  clutter 
problem  easier  than  that  scaled  factor  TFACF. 

The  output  SINR  of  full  adaptive  STAP  is  increased  by  8  dB  as  seen  in 
Fig.  4.12  (d).  This  increase  is  due  to  less  clutter  power  being  returned  which  in 
turn  decreases  the  main  beam  clutter  return.  The  small  main  beam  clutter  return 
increases  the  best  possible  SINR  close  to  the  thermal  noise  limit. 

4-5.2  Second  Gold  Code.  The  second  Gold  code  is  defined  as  [3] 

0101100101100001001000000111010.  (4.22) 

Again,  this  scheme  is  coded  by  using  a  binary  phase  scheme  within  the  chips.  A 
value  of  1  is  coded  as  e J7r  and  0  is  coded  as  0.  Figure  4.13  shows  the  TFACF  of  the 
second  Gold  code. 

Referencing  Fig  4.14  (a),  the  clutter  spectrum  of  the  gold  code  is  plotted. 
Examining  this  plot,  a  7  dB  degradation  of  the  clutter  return  power  is  gained  by 
implementing  this  compression  scheme.  This  gain  and  plot  is  the  same  as  the  first 
Gold  code.  Any  differences  are  minuet. 

The  minimum  variance  plot  of  the  second  gold  code,  Fig.  4.14  (b),  indicates 
that  resolution  capabilities  are  lost  by  using  a  compressed  pulse  width.  Comparing 
this  plot  against  Fig.  4.12  (b),  the  resolution  capability  is  very  similar  between  each 
Gold  code. 

Examining  the  bandwidth  of  the  second  Gold  code, 

31 

Bc  =  —  =  38.75  MHZ,  (4.23) 
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Figure  4.13:  The  Time  Frequency  Auto  Correlation  Function 
Plot  of  the  Second  Gold  Coded  pulse. 

an  increase  from  the  uncompress  pulse  bandwidth  by  more  than  37  MHZ  is  found. 
This  bandwidth  is  identical  to  the  first  Gold  code.  However,  equal  bandwidths  are 
expected  since  each  code  is  a  length  of  31.  Therefore,  the  thermal  noise  introduced 
into  the  system  is  identical  to  that  of  the  first  Gold  code.  The  eigenvalue  decompo¬ 
sition  shows  this  increase  in  thermal  noise,  Fig.  4.14  (c). 

The  output  SINR  of  full  adaptive  STAP  is  increased  by  8  dB  as  seen  in 
Fig.  4.14  (d).  This  increase  is  due  to  less  clutter  power  being  returned  which  in 
turn  decreases  the  main  beam  clutter  return.  The  small  main  beam  clutter  return 
increases  the  best  possible  SINR  close  to  the  thermal  noise  limit. 

4-6  Linear  Frequency  Modulation  Code 

The  last  waveform  examined  is  the  Linear  Frequency  Modulation  (LFM)  code. 
This  is  the  most  common  transmitted  pulse  and  is  frequently  used  with  Synthetic 
Aperture  Radar  systems.  As  outlined  in  Chapter  3,  the  LFM  is  not  discrete  like 
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(a)  (b) 


(c)  (d) 


Figure  4.14:  Second  Gold  Coded  Pulse 

(a)  Signal  match  describing  the  clutter  spectrum. 

(b)  MVE  describing  the  resolution  capability  utilizing  this  waveform. 

(c)  Eigenvalue  decomposition  illustrating  the  rank  of  the  clutter  spectrum  and  the 
overall  difficulty  of  the  problem. 

(d)  Output  SINR  of  fully  adaptive  Space  Time  Adaptive  Processing  (STAP)  illus¬ 
trating  target  detection  capability  at  all  normalized  Doppler  frequencies. 


Figure  4.15:  The  Time  Frequency  Auto  Correlation  Function 
Plot  of  the  Linear  Frequency  Modulated  pulse. 


the  other  waveforms  examined.  Therefore,  the  special  case  for  the  LFM  derived  in 
Chapter  3  is  used  to  model  the  effects.  To  ensure  that  the  derived  TFACF  in  chapter 
is  correct,  the  TFACF  of  the  LFM  is  required.  Figure  4.15  shows  the  results  from 
the  derived  TFACF.  This  figure  is  what  is  expected  for  the  TFACF  give  confidence 
that  the  derivation  is  correct. 

The  first  case  examined  for  the  LFM  case  is  the  signal  match  technique  as 
seen  in  Fig.  4.16  (a).  As  shown,  a  4  dB  loss  in  clutter  return  power  results  from 
implementing  the  LFM.  This  clutter  spectrum  resembles  the  other  clutter  spectrum 
found  using  the  discrete  waveforms. 

The  next  step  is  to  evaluate  how  the  resolution  is  effected  using  the  LFM. 
Figure  4.16  (b)  shows  that  range  resolution  is  lost  by  implementing  this  compressed 
waveform.  Any  targets  that  fall  into  the  main  beam  maybe  lost  due  to  high  clutter 
power. 
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(a)  (b) 


(c)  (d) 


Figure  4.16:  Linear  Frequency  Modulated  Coded  Pulse 

(a)  Signal  match  describing  the  clutter  spectrum. 

(b)  MVE  describing  the  resolution  capability  utilizing  this  waveform. 

(c)  Eigenvalue  decomposition  illustrating  the  rank  of  the  clutter  spectrum  and  the 
overall  difficulty  of  the  problem. 

(d)  Output  SINR  of  fully  adaptive  Space  Time  Adaptive  Processing  (STAP)  illus¬ 
trating  target  detection  capability  at  all  normalized  Doppler  frequencies. 
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The  bandwidth  for  the  LFM  waveform  is  calculated  differently  than  that  of 
the  discrete  waveforms.  Bandwidth  for  a  LFM  waveform  is  defined  as 

Bifm  =  h~h  (4.24) 

where  f\  and  /2  is  the  starting  and  stopping  frequency  from  the  chirped  waveform  [5] . 
For  this  simulation,  a  bandwidth  of  20  MHZ  was  used.  Since  this  waveform  is  smaller 
than  some  of  the  other  compression  schemes,  the  thermal  noise  introduced  to  the 
system  is  smaller.  Figure  4.16  (c)  shows  that  the  noise  floor  is  at  —131  dB.  This 
level  allows  more  clutter  returns  than  the  larger  compression  schemes  making  for  a 
difficult  clutter  problem. 

The  last  step  is  to  examine  the  output  SINR  of  the  fully  adaptive  STAP.  This 
output  is  shown  in  Fig.  4.16  (d).  A  gain  of  7  dB  is  achieved  from  implementing  this 
LFM  waveform.  This  gain  means  that  the  main  beam  clutter  is  not  as  large  as  it  is 
with  the  approximation.  Therefore,  targets  are  able  to  be  detected  with  more  ease. 

The  examination  of  several  different  waveforms  is  now  complete.  Analysis  was 
conducted  using  four  different  techniques,  signal  match,  minimum  variance  estima¬ 
tion,  eigenvalue  decomposition  and  output  SINR  of  fully  adaptive  STAP. 
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V.  Conclusions 


The  radar  model  encompasses  many  different  aspects  to  create  an  accurate  rep¬ 
resentation  of  real-world  environments.  To  complete  this  model,  an  under¬ 
standing  of  how  the  signal  is  transmitted  and  received  is  paramount.  The  received 
signal  contains  target  information.  In  a  perfect  world,  the  power  that  was  transmit¬ 
ted  would  be  the  same  power  received.  However,  there  are  outside  factors  such  as 
thermal  noise  and  ground  clutter  that  effect  radar  detection  capability  and  need  to 
be  accounted  for  in  the  model. 

Typically,  the  radar  model  makes  an  approximation  on  processing  the  received 
data.  This  approximation  is  that  there  is  no  Doppler  or  range  shift  when  matched 
filtering  the  received  signal  with  the  transmitted  signal.  This  assumption  makes 
the  Time  Frequency  Auto  Correlation  Function  (TFACF)  equal  one.  This  research 
removes  that  assumption  and  analyzes  the  performance  effects  on  the  radar  model 
using  the  TFACF. 

After  deriving  the  generic  TFACF  expression,  simulations  were  run  for  several 
different  waveforms  including  a  13-length  Barker,  26-length  concatenated  Barker,  two 
Frank  codes  (M=5  and  M=10),  two  Gold  codes  and  a  Linear  Frequency  Modulated 
code.  From  the  results  seen  in  Chapter  4,  there  were  some  performance  degradations 
and  gains  on  the  radar  model.  The  performance  degradation  was  apparent  within 
the  resolution  capabilities.  Each  waveform  experienced  a  loss  in  resolution  which 
overall  effects  the  radar’s  ability  to  distinguish  targets  located  close  to  each  other. 

Improvements  were  noticed  within  the  reduced  amount  of  returned  clutter 
power  as  well  as  the  difficulty  of  the  clutter  problem.  All  the  implemented  waveforms 
reduced  the  power  return  from  the  clutter.  This  reduction  allows  for  better  target 
detection  since  the  target  power  has  less  interference  to  compete  against.  Thermal 
noise  plays  a  major  role  with  in  this  reduction.  The  amount  of  thermal  noise  allowed 
is  directly  proportional  to  the  bandwidth  of  the  pulse.  The  larger  the  bandwidth, 
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the  greater  effect  from  thermal  noise.  This  increase  in  thermal  noise  dominates 
the  clutter  power  being  returned.  Thus,  the  clutter  power  is  not  as  effective  on 
the  system.  This  is  also  seen  in  the  eigenvalue  decompositions.  The  Eigenvalue 
decomposition  provides  a  metric  describing  the  difficulty  of  the  clutter  environment. 
In  all  the  decomposition  plots,  the  thermal  noise  floor  is  raised  considerably.  This 
raised  effect  dominates  less  clutter  return  values  making  the  problem  less  difficult. 

Like  every  radar  problem,  there  are  trade-offs.  By  approximating  the  TFACF 
to  one,  there  are  performance  characteristics  that  are  not  accounted  for  that  effect 
the  system.  Assuming  this  approximation,  the  radar  model  provides  better  range 
resolution  than  what  is  actually  achievable.  However,  approximating  the  TFACF,  a 
worse  case  scenario  is  established.  This  scenario  allows  for  a  more  difficult  clutter 
power  and  more  main  beam  clutter  power  introduced  into  the  model.  This  intro¬ 
duction  hampers  target  detection  capabilities  with  the  radar  model.  Therefore,  a 
decision  must  be  made  when  creating  the  radar  model.  If  accurate  resolution  is 
required,  the  TFACF  approximation  is  not  recommended.  However,  if  the  model  is 
created  to  overcome  clutter  power,  then  the  approximation  is  fine  due  to  the  over 
estimation  of  actual  clutter  power. 

By  assuming  the  TFACF  is  approximately  one,  important  information  is  lost. 
Depending  on  the  waveform  and  its  bandwidth,  the  picture  of  the  interference  envi¬ 
ronment  could  be  distorted.  There  is  some  error  within  the  current  assumption  but 
not  enough  to  include  in  the  radar  model. 

Future  work  with  this  topic  could  include  removing  the  limitation  of  the  pulse 
only  mismatching  by  less  than  one  half  a  chip  or  how  multiple  time  around  pulses 
could  effect  the  model.  Another  idea  is  to  integrate  a  barrage  noise  jammer  into  the 
model  to  see  the  effects,  if  any,  the  jammer  has  on  this  model. 
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Appendix  A.  Matlab  Code 


This  appendix  contains  the  Matlab®  code  used  to  model  both  the  conventional 
and  the  scaled  Time  Frequency  Auto  Correlation  Function  (TFACF)  radar 
models. 


A.l  TFACF  Function  (Discrete  Case) 

Listing  A.l:  A  function  calculating  the  discrete  case  of  the  TFACF. 

(appendixl/tfacf.m) 

%============================================================== 

7.  Created  By:  Capt  Gregg  Izdepski 
7.  Date:  21  Mar  05 
7. 

"/.This  function  calculates  the  discrete  form  of  the  Time  Frequency 
‘/.Auto  Correlation  Function.  There  is  the  capability  to  implement 
'/.seven  different  pulse  compression  waveforms  within  the  code. 
“/.Simply  remove  the  comments  of  the  desired  code  and  run  the 
“/.program.  The  variable  ‘waveform1  is  passed  for  the  file  name 
“/.while  saving  the  plots.  ‘Tc’  is  passed  to  help  calculate  the 
"/.bandwidth,  ‘chips’  contains  the  value  of  the  TFACF. 
%============================================================== 


function  [chips ,Tc ,waveform]=TFACF (tau ,f_t , numSubCell s , q) 
clear  chip ; 

7.  %====================================================== 

7,  7.  simple  rect  pulse 
7.  psi  =  zeros  (1,1)  ; 

7,  wavef  orm=  ’  rect  ’  ; 


7.  %========================================== 

7,  "/.13-length  Barker  Code 

psi= [zeros (5 , 1) ;  pi;  pi;  0;  0;  pi;  0;  pi;  0] ; 

waveform= ’ Barker ’ ; 

7,  psi  =  [psi  ;  f  lipud  (psi)  ]  ; 

7,  wavef  orm=  ’  Con.Barker  ’  ; 


7.  1  =  =  =  =  =  =  =  =  =  = 
7.  “/.Frank  Code 
7.  M  =  10 ; 

7. 
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7. 
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7. 

7. 

7. 
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7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 


7.  This  is  for  values  in  degrees 
x=linspace (0 , 360-360/M, M)  ; 
y=  [0:1: M-l]  1 ; 
matr ix=y*x ; 

max_value=round ( (max (max (matrix) ) ) /360) ; 

7.  This  is  for  values  in  radians 

x=linspace  (0 , 2*pi -2*pi/M  ,  M)  ;  7»for  values  in  radians 
y=  [0  :  1  :  M-l]  »  ; 
matr ix=y*x ; 

max_value=round ( (max (max (matrix) ) ) /2/ pi) ; 

7.  Creating  the  matrix 
for  p=l : max.value 
for  k=l : M 

temp_matrix= [matrix (k , :)] ; 
for  m=l : M 

7. This  is  for  values  in  degrees 

if  temp.matrix ( : ,m)>=360 

temp.matrix ( : , m) =temp_matr ix ( : , m) -360 ; 

else 

temp_ matrix ( : , m) =temp_matr ix ( : , m) ; 

end 

7 This  is  for  values  in  radians 

if  temp.matrix ( : ,m)>=2*pi 

temp.matrix ( : , m) =temp_matr ix (:,m)-2*pi; 

else 

temp.matrix ( : , m) =temp_matr ix ( : , m) ; 

end 

end 

interim(k, :) =temp_matr ix ; 

end 

matrix=interim ; 

end 


i  psi=matrix ( : ) ; 

7.  waveform=’frank_10  1  ;  '/.check  that  the  #  equals  M 


/.  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  = 
'/.Gold  Code 

7.  code_l  =  [0001101111011010001111110100000]  ’*pi; 

1  code_2=  [0101100101100001001000000111010]  ’*pi; 

1 

%  psi  =  code_l;  '/.change  the  _#  to  select  code  #1  or  #2 

i  wavef  orm=  ’  gold_  1  '/.check  that  numbers  are  matching 
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P=length (psi ) ; 
Tc=tau/P ; 


7,  #  of  chips 
%  Chip  time 


/.  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  = 
90  7«subdivsions  within  the  ring 

’/  matching  the  document  ...  tau.prime  is  the  delay,  tau-T_R 
tau_prime=linspace (-Tc/2 , Tc/2 , numSubCells ) ; 
tau_prime=tau_prime (q) ; 

95 

7.  =  =======  =  =  =  =  =  =  =  ====  =  =  =  =  =  =  =======  =  =  =  =  =  =  =  ===  =  =  =  =  =  =  =  ======= 

7o  calculat  ing  the  TFACF 

chips = zeros ( length (tau.prime) ,1) ; 

100  7.  The  less  than  case  expression 
if  tau.prime  <  0 
for  p=0 : P-1 

if  p~=P-l 

chip (p+1) =1/P/Tc * (  ... 

105  (Tc+tau_pr ime ) * . . . 

exp(j*pi*f_t*((2*p+l)*Tc+tau_prime))*. . . 
sinc(f_t*(Tc+tau_prime))-. . . 
exp(j*psi(p+l)-j*psi(p+2))*tau_prime*. . . 
exp(j*pi*f_t*(2*(p+l)*Tc+tau_prime))*. . . 
110  s inc ( f _t * t au_pr ime ) )  ; 

else 

chip (p+1) =1/P/Tc * (  ... 

(Tc+tau_pr ime ) * . . . 

exp(j*pi*f_t*(Tc*(2*p+l)+tau_prime))*. . . 
115  sine (f_t*(Tc+tau_prime) ) ) ; 

end 

end 

chips  =  sum ( sum ( chip )) ; 

120  clear  chip  ; 

7.  The  greater  than  case  expression 
elseif  tau.prime  >  0 
for  p=0 : P-1 
125  if  p==0 

chip (p  +  1) =1/P/Tc *  (  ... 

(Tc - 1 au_pr ime ) *  .  .  . 

exp(j*pi*f_t*((2*p  +  l)*Tc  +  tau_prime))*.  .  . 
sine (f_t*(Tc-tau_prime) ) )  ; 

130  else 

chip (p+1) =1/P/Tc * (  ... 

(Tc-tau_prime)*. . . 

exp(j*pi*f_t*((2*p+l)*Tc+tau_prime))*. . . 
s inc ( f _t * ( Tc - t au_pr ime ) ) + . . . 

135  exp(j*(psi(p+l)-psi(p)))*tau_prime*... 

exp ( j *pi*f _t * (2*p*Tc  +  tau.prime ) ) * . . . 
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end 


s inc ( f _t * t au_pr ime ) )  ; 


end 

140 

chips  =  sum ( sum ( chip ))  ; 
clear  chip ; 


else 

145 

'/.The  equal  to  case  expression 
chips  =  s inc ( f _t *p*Tc ) ; 

end 


A. 2  TFACF  Function  (LFM) 

Listing  A. 2:  A  function  calculating  the  continuous  LFM  case  of  the  TFACF. 

(appendixl /TfacfLfm.m) 


7.============================================================== 

7,  Created  By:  Capt  Gregg  Izdepski 
7.  Date:  21  Mar  05 
1 

5  7»This  function  calculates  the  Time  Frequency  Auto  Correlation 
7.Function  of  the  continuous  LFM  waveform.  It  returns  the  TFACF 
7»value  in  'chips’  and  ‘waveform’  is  the  root  filename  for  saving 
7.the  analysis  plots. 

/.  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  = 

10 

function  [chips ,b , waveform] =TFACF_lfm (tau , f_t , numSubCells , q) 

waveform=’lfm’ ; 

15  b  =  -5el3;  7.chirp  rate 

tau_prime=linspace (-tau/2 , tau/2 , numSubCells ) ; 

tau_prime=tau_prime (q) ; 

chips  =  zeros ( length (tau.prime)  ,1)  ; 

20 

7.  The  less  than  expression 
if  tau.prime  <  0 

chip=(l/tau)*exp(-j*2*pi*b*tau_prime~2)*. . . 

exp(j*pi*(2*b*tau_prime+f_t)*(tau+tau_prime))*. . . 

25  ( t au  + 1 au_pr ime ) * .  .  . 

sine ((2*b*tau_prime  +  f_t)*(tau  +  tau_prime))  ; 

chips = sum ( sum ( chip ) ) ; 
clear  chip ; 

30 

7.The  greater  than  expression 
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elseif  tau.prime  >  0 

chip=(l/tau)*exp(-j*2*pi*b*tau_prime~2)*.  .  . 

exp ( j *2* pi  * ( 2* b * t au_pr ime  +  f _t ) * t au_pr ime ) * .  .  . 

35  exp(j*pi*(2*b*tau_prime+f_t)*(tau-tau_prime))*. . . 

(tau-tau_prime) * . . . 

sine ((2*b*tau_prime+f_t)*(tau-tau_prime)) ; 

chips = sum ( sum ( chip ) ) ; 

40  clear  chip  ; 

'/.The  equal  to  expression 
else 

chips  =  sine ( f _t *tau) ; 

45  end 


A. 3  Antenna  Pattern  Function 

Listing  A. 3:  A  function  calculating  the  necessary  antenna  pattern  for  weighting 

the  clutter  returns. 

(appendixl  /  antenna,  m) 


%=========================================================== 

•/.  Created  by:  Capt  Gregg  Izdepski 
/  Date:  21  Mar  05 

7. 

5  7,  This  function  builds  the  transmitted  antenna  pattern.  The 
7.  pattern  built  here  is  used  to  scale  the  clutter  returns. 
%=========================================================== 


function  [g ,  g_log ,  w_af ] =antenna (theta ,  phi,  N,  P,  fo) 

10 

7oGlobal  Parameters 

c=3e8;  7»  speed  of  propogation 

lambda  =  c/fo;  7»  wavelength 

15  7»backlobe  attenuation  factor 

dB=-30;  /  dB  attenuation  factor 

back  =  10~  (  dB  /  1 0 )  ;  7«  backlobe  weighting  factor 

7«  building  the  element  pattern 
20  g=(abs(cos(theta)*cos(phi’))) .“2; 

a=[find(phi  <=  -90*pi/180);  find(phi  >=  90*pi/180)]; 
temp=g(: ,a)*back; 
g ( :  , a) =temp  ; 

25  g_log=10* loglO ( abs (g/max (max (g) ) ) ) ;  7.normalizing  g 

70raising  the  floor  to  -50  dB 

g_log  =  (g_log  >  -50)  . *  g_log  +  ( g  _ 1 o  g  <=  -50) * ( -50)  ; 
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30  '/.array  factor 

delta_z=lambda/2 ; 
delta_x=lambda/2 ; 

var  _x  =  (  pi  *  delt  a_x  /  lambda)  ;  '/,dx  portion  of  array  factor 
var  _z  =  (  pi  *  delt  a_z  /  lambda)  ;  °/.dz  portion  of  array  factor 
35 

7.  Building  the  Array  Factor 
for  i=l : length (theta) 

for  h=l : length (phi ) 

w_af(i,h)=exp(j*var_z*(P-l) . *sin(theta(i))) .  *  . . . 

40  exp(j*var_x*(N-l).*cos(theta(i)).*sin(phi(h)))... 

,*sin(var_z*P.*sin(theta(i))) . *sin(var_x*N.*. . . 
cos(theta(i))  . *  s in ( phi (h) ) )  . / .  .  . 
(sin(var_z.*sin(theta(i))) .*sin(var_x.*. . . 
cos(theta(i))  . *  s in ( phi (h) ) ) )  ; 

45  end 

end 

7,  Array  Factor  special  cases 
w_zeros=N*P ; 

50  w_thet a  =  N* exp ( j * var _z * (P - 1 )  . * s in ( thet a ) )  .  . 

((sin(var_z*P.*sin(theta))) ,/(sin(var_z.*sin(theta)))) ; 
w_phi=P*exp(j*var_x*(N-l) .*sin(phi ’ ) ) • • 

((sin(var_x*N.*sin(phi’)) ,/(sin(var_x.*sin( phi’ ))))); 

55  7, the  complete  array  factor  matrix  (accounting  for  special  cases) 
if  ~ i sempty ( f ind ( phi ==0) ) 

w_af ( : , f ind  (phi ==0) )=w_theta; 

end 

60  if  ~ i sempty ( f ind ( theta ==0) ) 

w_af ( f ind ( thet a  =  =  0 )  ,:)=w_phi; 

end 

if  ~isempty  ( [f ind (theta==0)  f ind (phi ==0) ] ) 

65  w_af(find(theta==0),f ind ( phi ==0) )=w_zeros ; 

end 

return 


A. 4  Conventional  Clutter  Model  Function 

Listing  A. 4:  A  function  calculating  the  conventional  clutter  model, 

(appendixl  /  clutter .  m) 


%================================ 

7«  Created  by:  Capt  Gregg  Izdepski 
7.  Date:  21  Mar  05 
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7. 


7«This  function  calculates  the  clutter  covariance  matrix  for  the 
7»  convent  ional  model.  Within  this  model,  the  Time  Frequency  Auto 
/•Correlation  Function  is  assumed  to  equal  one. 


7. 


function  [R_clut]=clutter(M,N,P,fo ,prf ,alt ,tau ,R,B,delta_x . . . 
,delta_z ,Pt ,elem_gain ,Ls) 


/•  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  = 
/•Known  parameters 
r_earth=63Tle3 ; 
c=3  e8  ; 

gamma_dB  =  -3 ; 
Nc=361 ; 

To=290 ; 

k=l . 380685e-23 ; 


7o radius  of  the  earth  in  m 
7o prop ogat  ion  speed 
/•clutter  gamma  in  dB 

7, #  of  clutter  patches  around  the  clutter  ring 
7,  st  andard  temp  (  K) 

7,  Boltzman  constant 


/•  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  = 
/•Calculated  parameters 
Tr=l/prf ; 

K=f loor (Tr/tau) ; 
a_earth=4*r_earth/3; 
gamma= 10 ~ ( gamma_dB / 1 0 ) ; 
velocity=delta_x/ (2*Tr ) ; 
lambda=c/ f o ; 

No=k*To ; 
sigma=B*No ; 


7.PRI 

7,  number  of  range  cells 
7,  effective  radius  of  the  earth 
7. gamma  in  watts 
7,  a  /  c  velocity  (m/s) 

7» lambda  in  m 
7. standard  noise  No 
7»sigma  for  the  noise 


/•  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  = 
70  Clutter  Covariance  Matrix 

Ru  =  c/  (prf  *2)  ;  7«unambigous  range 

Rh=  sqrt  (  alt  ~  2  +  2*  alt  *  a.earth  )  ;  70horizontal  range 

7«determining  Rc  range  for  each  bin 
bin_time=(Tr-tau) /K ; 

range=(0: c*bin_time/2: c*K*bin_time/2-c*bin_time/2) ; 

70locating  the  range  bin  for  the  range 
bin_desired=find (abs (range -R)==min(abs (range -R) ) ) ; 

Rclut  =  range  (bin_desired)  ;  7«cycling  current  Re's 

Nr  =  Rh/Ru;  7,  calculat  ing  Nr 

if  Rclut  *Nr  <=  Rh 
Nr  =  ceil  (Nr )  ; 

else 

Nr=f loor  (Nr) ; 

end 
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Ri  =  (0  :  Nr  -  1 )  * Ru  +  Rclut  ; 
del_r=c/ (2*B) ; 


Zambigous  range  cells 
‘/.range  resolution 


'/.azimuth  clutter  angles 
del_phi=2*pi/Nc ; 

phi_clut=(-pi  :  del_phi  :pi-del_phi)  ’  ; 

%  elevation  clutter  angles 

theta_clut=-asin(((Ri.~2)+alt*(alt+2*a_earth)) ./. . . 
(2*Ri.*(alt+a_earth))) ’ ; 

'/.creating  the  antenna  pattern  (antenna  fn) 

[g  ,  g_log ,  w_af ]  =  antenna (theta_clut  ,  phi_clut  ,  N,  P,  fo); 

% grazing  angle 

psi_clut=-asin(((Ri.~2)-alt*(alt+2*a_earth)) ,/(2*Ri.*a_earth)); 

% int ial iz ing  the  covariance  matrix 
R_c=zeros (M*N*P ,M*N*P) ; 

vis  =  waitbar(0,  'Calculating  Rings  ...  1 )  ; 
for  i=l : Nr 

if  Ri ( i ) >alt ; 


7»  ground  reflectivity 

s igma_ c =gamma* s in ( ps i_c lut ( i ) ) .*Ri(i) .*del_phi*. . . 
del_r . *  sec (psi_clut (i) ) ; 

for  k=l : length (phi_clut ) 

var _x  =  ( delt a_x * s in ( phi _clut (k) ) *  .  .  . 

cos (theta_clut (i) ) ) / lambda ; 
var_z=(delta_z*sin (theta.clut (i) ) ) / lambda ; 
w_bar=(2*velocity*cos(theta_clut(i))*. . . 
sin(phi_clut (k) )/ lambda) *Tr ; 

steer_a=exp ( j *2*pi *  [0 : N-l]  ’*var.x)  ; 
steer_b=exp ( j *2*pi * [0 : M-l]  ’*w_bar)  ; 
steer_e=exp ( j *2*pi * [0 : P-1]  ’*var.z)  ; 

7  the  steering  vector 

v=kron(steer_e ,kron(steer_b , steer.a) ) ; 

G=abs(w_af (i,k) . "2)*g(i,k) ; 

7.CNR  per  patch 

xsi_c (i ,k)  =  (Pt*G*g(i ,k)*elem_gain*( lambda  ~2) * .  .  . 
sigma_c) / ( ( (4*pi) ~3) *No*B*Ls*Ri (i) "4) ; 

7building  the  covariance  matrix 
R_c=R_c+sigma*xsi_c (i,k)*v*v'; 
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end 


end 

waitbar (i/Nr , vis) 

end 

110  close  (vis) 

R_  clut  =  R_c ; 


A. 5  Scaled  Clutter  Model  Function 

Listing  A. 5:  A  function  calculating  the  TFACF  scaled  clutter  model, 
(appendixl  /  ClutterDop.m) 


%============================================================== 

"/.  Created  by:  Capt  Gregg  Izdepski 
"/.  Date:  21  Mar  05 

7. 

5  "/.This  function  calculates  the  clutter  covariance  matrix  for  the 
70model  scaled  by  the  Time  Frequency  Auto  Correlation  Function 
"/.(TFACF)  .  The  code  in  this  model  is  very  similar  to  that  of  the 
"/.conventional  clutter  model.  The  difference  is  noted  when 
"/.finding  the  steering  vector  b.  Steering  vector  b  is  scaled  by 
10  "/.the  appropriate  TFACF  value. 

»/„  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  = 
function  [R_clut_dop,  sigma,  wavef orm]  =  clut t er _dop (M , N , P , f o , prf . .  . 

f  •  •  • 

alt ,tau,R,B,delta_x ,delta_z ,Pt , elem.gain ,Ls , numSubCells ) 


15  "/.  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  = 
"/.Known  parameters 
r_earth=6371e3 ; 
c=3  e8 ; 

gamma_dB  =  -3 ; 

20  Nc  =36 1  ; 

To=290  ; 

k=l . 380685e-23 ; 


7, radius  of  the  earth  in  m 

“/.progation  speed 
7, clutter  gamma  in  dB 

°/0#  of  clutter  patches  around  the  clutter  ring 
7, standard  temp  (K) 

'/.Boltzman  constant 


25  °/0  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  = 
'/.Calculate  parameters 
Tr=l/ prf ; 

K=f loor (Tr/tau) ; 
a_earth=4*r_earth/3; 

30  gamma= 10 ~ ( gamma_dB / 1 0 ) ; 
velocity=delta_x/ (2*Tr ) ; 
lambda  =  c/ f o  ; 

No=k*To ; 
sigma=B*No ; 

35 

/.  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  = 


y.pRi 

'/.number  of  range  cells 
°/0  effective  radius  of  the  earth 
’/.gamma  in  watts 
°/0a/c  velocity  (m/s) 

“/.lambda  in  m 
"/.standard  noise  No 
"/.sigma  for  the  noise 
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7.  Clutter  Covariance  Matrix 


Ru=c/ (prf  *2)  ;  '/.unambigous  range 

Rh  =  sqrt  (  alt  ~  2  +  2*  alt  *  a.earth  )  ;  ‘/.horizontal  range 

'/.determining  Rc  range  for  each  bin 
bin_time=(Tr-tau) /K ; 

range=(0: c*bin_time/2: c*K*bin_time/2-c*bin_time/2) ; 

'/.locating  the  range  bin  for  the  range 
bin_desired=find (abs (range -R)==min(abs (range -R) ) ) ; 

Rclut  =  range  (bin_desired)  ;  “/.cycling  current  Re’s 

Nr=Rh/Ru;  “/.calculating  Nr 

7.  Nr  =  1 ; 

if  Rclut  *Nr  <=  Rh 
Nr  =  ceil  (Nr )  ; 

else 

Nr=f loor  (Nr) ; 

end 


Ri  =  (0  :  Nr  -  1 )  *  Ru  +  Rclut;  '/.ambigous  range  cells 

del_r=c/ (2*B) ;  '/.range  resolution 

"/.azimuth  clutter  angles 
del_phi=2*pi/Nc ; 

phi_clut=(-pi : del_phi :pi-del_phi) 1 ; 

7.  int  ial  iz  ing  the  covariance  matrix 
R_c=zeros(M*N*P,M*N*P) ; 

"/.Dividing  each  main  range  cell  into  subcells 
r_sub=linspace (0 , del_r , numSubCells ) ; 

vis  =  waitbar (0 ,  “Calculating  Rings  ...  ’)  ; 
for  i = 1 : Nr 

if  Ri ( i ) >alt ; 

disp([“Range  ring  center  is  ’  ,  num2st r ( Ri ( i ) / 1 e3 )  ,  ’  km’]); 

vis_a  =  waitbar  (0  ,  1  Calculating  Sub rings  .  .  .  ’  )  ; 
for  q  =  0 : numSubCells -1 

“/.fining  range  to  each  smaller  cell 
Rq=Ri ( i ) -del_r /2+ r _sub ( q+ 1 ) ; 

7,  grazing  angle 

ps  i  _clut  =  - as  in  (  (  (  Rq  “  2  )  -  alt  *  (  alt  +2*  a.earth  )  )  /  .  .  . 
(2*Rq*a_earth) )  ; 
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°/0  ground  reflectivity 

s igma_c lut =gamma* s in ( ps i_c lut ) *Rq* del_phi * . . . 
mean(diff (r_sub))*sec(psi_clut) ; 

'/.elevation  clutter  angles 

thet  a_c  lut  =  - as  in  (  (  (  Rq  ~  2)  +  alt  *  (  alt  +  2*  a.earth  )  )  /  .  .  . 
(2*Rq* ( alt+a_earth) ) )  ; 

var_z=(delta_z*sin (theta.clut ) ) / lambda ; 
steer_e  =  exp ( j *2*pi *  [0 : P-1]  ’*var.z)  ; 

for  k=l : length (phi_clut ) 

'/.creating  the  antenna  pattern  (antenna  fn) 
[g,g_log ,w_af]  =  antenna (theta_clut  ,phi_clut(k)  ,  .  .  . 
N  ,  P  ,  f  o  )  ; 

var_x=(delta_x*sin(phi_clut(k))*. . . 

cos(theta_clut))/ lambda ; 
w_bar = (2* ve loc i ty *  cos ( thet a_ clut ) * .  .  . 
sin (phi_clut (k) ) / lambda ) *Tr ; 

'/.These  next  two  lines  determine  which  TFACF  is  called.  One  of 
110  '/.these  two  should  be  commented  out  while  running  sims  . 

'/.Otherwise  ,  only  the  LFM  TFACF  will  be  return  (since  it  is  last) 

7,  This  is  for  the  discrete  TFACF 

[chips  ,Tc,waveform]=TFACF(tau,w_bar/Tr,  .  .  . 

115  numSubCells , (q+1) ) ; 

B  =  l/Tc;  /.calculates  the  bandwidth 

'/.This  is  for  the  LFM  TFACF 

/  [chips,  b,  wavef orm] =Tf acf Lfm (tau , w_bar/Tr , . . . 

120  /  numSubCells , (q+1) ) ; 

l  B  =  b*tau;  7.LFM  bandwidth 

sigma=B*No;  “/,  calculat  ing  the  noise  power 

125  steer_a=exp(j*2*pi*[0:N-l] ’ *var_x) ; 

steer_b=exp(j*2*pi*  [0 : M- 1]  ’*w_bar)*chips ; 

v_dop=kron(steer_e , (kron(steer_b , steer.a) ) ) ; 

130  G=abs (w_af ~2) *g ; 

xsi_c  =  (Pt*G*g*elem_gain*( lambda  ~2) *  s igma_clut ) . .  . 

/  .  .  . 

( ( (4* pi ) ~3) *No*B*Ls  *Rq ~4)  ; 

135  R_c=R_c+sigma*xsi_c*v_dop*v_dop ’ ; 

end 
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waitbar ((q+1) / numSubCells , vis_a) ; 


end 

140  close  (vis_a); 

end 

waitbar(i/Nr,vis) ; 


end 

145  close  (vis); 

R_clut_dop=R_c ; 


A. 6  Overall  Model 

Listing  A. 6:  The  overall  model  incorporating  all  functions.  Analysis  is  completed 

here  allowing  for  analysis. 

(appendixl /ModelClutter.m) 

7.============================================================== 

%  Created  by:  Capt  Gregg  Izdepski 
7.  Date:  21  Mar  05 
% 

5  7oThis  code  incorporates  all  f  unct  ions  to  create  the  two  different 
7»radar  models,  conventional  and  scaled.  This  code  should  remain 
7»unchanged  since  all  parameters  are  changed  within  the  functions. 

7»The  only  changes  are  title  descriptions,  if  desired.  The  plots 
7»are  directly  saved  to  the  current  file  directory 
10  7»(both  *.*eps  and  *.*png). 

*/.  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  = 
close  all;  clear  all;  clc ;  format  compact; 

tic; 

15  /.=  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  = 

7oGlobal  Parameters 
M  =  16  ; 

N  =  16  ; 

P  =  2; 

20  f o  =  l . 24 e9 ; 
prf =1984; 
t  au  =0 . 8 e -6  ; 

Pt  =  1 . 5  e3  ; 

B=l/ tau ; 

25  Fn_dB  =3 ; 
alt  =9  e3  ; 

R  =  66  e3 ; 

delta_x=0 . 10922; 
delta_z  =0 . 140716 ; 

30  Ls_dB  =3 ; 

elem_gain_dB=0; 

To=290 ; 

k=l . 380685e -23 ; 


70  pulses 

7o#  of  azimuth  channels 
7«#  of  elevation  channels 
7.  transmitted  frequency  (GHz) 
7«pulse  repetition  frequency 
7tpulse  width 
7o  power  transmitted  (kW) 

7«bw  in  Hz 

7.noise  figure  in  dB 
70a/c  altitude  (m) 

7«Range  to  target  in  m 
7«azimuth  channel  spacing  (m) 

7«  elevation  channel  spacing  (m) 
7«system  loss  in  dB 
7« element  gain  in  dB 
7«  st  andard  temp  (  K) 

7oBoltzman  constant 
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c=3e8;  '/.propogation  speed 

velocity=90;  °/«a/c  velocity  (m/s) 

numSubCells  =  101  ;  "/#  of  smaller  clutter  rings  in  each  main  ring 


"/  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  = 
"/calculated  constants 

Ls  =  1 0  ~  (  Ls_dB  /  10  )  ;  "/system  loss  in  watts 

elem_gain  =  1 0 “  (  elem_gain_dB  /  10 )  ;  "/element  gain  in  watts 
No=k*To;  "/standard  noise  No 

sigma  =  B*No;  "/sigma  for  the  noise 

lambda  =  c/fo;  "/lambda  (m) 


"/============================================================== 

"/Clutter  Covariance  Matricies 

“/Clutter  assuming  no  doppler  tolerance 

[R_clut]=clutter(M,N,P,fo,prf  ,alt  , tau,R,B,delta_x ,delta_z  ,Pt  ,  .  .  . 
elem_gain , Ls ) ; 

"/Clutter  w/  doppler  effect 

[R_clut_dop ,sigma_dop ,waveform]=ClutterDop(M,N,P,fo,prf ,alt , . . . 
tau ,R,B,delta_x ,delta_z ,Pt ,elem_gain ,Ls , numSubCells ) ; 


"/  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  = 

"/Thermal  Noise  Block  Matrix 

R_noise=sigma*eye(M*N*P) ; 

"/calculating  new  noise  matrix  for  implement  waveform 
sigma_dop=B_dop*No ; 

R_noise_dop=sigma_dop*eye (M*N*P) ; 


"/  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  = 
"/Total  R  matrix  (including  clutter  and  thermal  noise) 

R_t ot  =  R_ clut +R_noi se  ;  "/conventional  data 

R_tot_dop  =  R_clut_dop  +  R_noise_dop  ;  "/scaled  data 


"/  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  = 
"/creating  the  steering  vector 

var_x=(delta_x*sin (0) *cos (0) ) / lambda ; 
var_z=(delta_z*sin (0) ) /lambda ; 
pad=2~ ( log2 (M) +2) ; 

w_bar  =  linspace (-0.5 ,0.5-1/pad ,pad)  ’ ; 

steer_a=exp(j*2*pi*[0:N-l] ’*var.x) ; 
steer_b=exp(j*2*pi*[0:M-l] 1 *w_bar ’) ; 
steer_e=exp(j*2*pi*[0:P-l] ’ *  var _z ) ; 


v=kron(steer_e ,kron(steer_b , steer.a) ) ; 
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%============================================================== 

"/.Spectra  Analysis 

%  S ignal  Match  ( SM) 

P_sm=diag(v’*R_tot*v)  ./diag(v’*v*R_tot(l,l))  ; 
P_sm_dop=diag(v’*R_tot_dop*v) ,/diag(v’*v*R_tot_dop(l,l)) ; 

f igur e  ; 

plot  (w_bar , 10*logl0 (abs (P_sm) ) ,w_bar , 10*logl0 (abs (P_sm_dop) ) . . . 

, ’ r-- ’ ) ; grid  on ; 

xlabel ( ’ Normalized  Doppler  (Hz)’); 
ylabel (’ Magnitude  (dBw)’); 

legend  ( 1 AF  Approximation ’,’ AF  Scale  Factor’); 

7o  title  (’Signal  Match  (SM)  Method’); 

"/„  saving  the  SM  plotted  results 
print(’-depsc’ .[waveform  ’_sm_dop.eps’]) ; 
print(’-dpng’ , [waveform  ’_sm_dop.png’]) ; 

'/.Minimum  Variance  Estimator  (MVE) 

P_mve  =  l./diag(v’*(R_tot\v))  ; 

P_mve_dop=l . / diag (v ’ * ( R_tot_dop\v) ) ; 

f igur e  ; 

plot (w_bar , 10*logl0 (abs ( P_mve ) /max ( abs ( P_mve ))) ,w_bar , . . . 

10*logl0 (abs (P_mve_dop) /max (abs ( P_mve_dop ))),’r--’);grid  on; 
xlabel (’ Normalized  Doppler  (Hz)’); 
ylabel (’ Magnitude  (dB)’); 

legend  ( ’ AF  Approximation ’,’ AF  Scale  Factor’); 

"/„  title  (’Minimum  Variance  Estimator  (MVE)  Method’); 

7  saving  the  MVE  plotted  results 
print ( ’ -depsc ’ .[waveform  ’_mv.eps’]) ; 
print ( ’ -dpng ’  ,  [waveform  ’ _mv . png  ’  ] )  ; 

'/.Eigenvalue 
Eig_R_tot=eig (R_tot ) ; 

Eig_R_t ot _dop = e ig ( R_t ot _dop ) ; 

f igur e  ; 

plot(l:length(Eig_R_tot)  ,  10*logl0( abs (sort (Eig_R_tot)))  ,  ... 

1 : length (Eig_R_tot_dop)  ,  10*logl0( abs (sort (Eig_R_tot_dop)))  .  .  . 
,  ’  r--  ’  )  ; grid  on ; 

xlabel (’ Sorted  Eigenvalue  Number’); 
ylabel (’ Eigenvalue  Magnitude  (dBw)’); 
legend(’AF  Approximation’.’AF  Scale  Factor ’,2); 
i  t itle (’ Eigenvalue  Comparison’); 
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7,  saving  the  Eigenvalue  plotted  results 
print (’-depsc’ , [waveform  ’_eig.eps’]) ; 
print  ( ’ -dpng ’  ,  [waveform  ’ _eig . png ’ ] )  ; 

140 

'/.Optimum  SINR 

SINR=diag ( s igma* v ’ * ( R_t ot \ v) ) ; 

SINR_dop=diag ( sigma_dop*v ’*(R_tot_dop\v)) ; 

145 

f igur e  ; 

plot(w_bar  ,  10*logl0(abs(SINR))  ,w_bar  , 10*logl0(abs(SINR_dop))  . 
, ’ r-- ’ ) ; grid  on ; 

xlabel ( ’ Normalized  Doppler  (Hz)’); 

150  ylabel (’ Magnitude  ( dB ) 1 ) ; 

legend/ ’AF  Approximation’,’AF  Scale  Factor’, 4); 

7,  t  itle  (’  Output  SINR  for  a  Fully  Adaptive  STAP  Process’); 

7.  saving  the  SINR  plotted  results 
155  pr int ( ’ -depsc ’ , [wavef orm  ’_sinr.eps’]); 
print/’- dpng’  ,  [waveform  ’_sinr. png  ’  ] )  ; 


time=toc / 60 
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